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Abstract 

We study the entropy of entanglement of the ground state in a wide family of one- 
dimensional quantum spin chains whose interaction is of finite range and translation 
invariant. Such systems can be thought of as generalizations of the XY model. The 
chain is divided in two parts: one containing the first consecutive L spins; the second 
the remaining ones. In this setting the entropy of entanglement is the von Neumann 
entropy of either part. At the core of our computation is the explicit evaluation of 
the leading order term as L ^ oo of the determinant of a block- Toeplitz matrix with 
symbol 

i\ g{z) 
g-\z) iX 



$(z) 



where g{z) is the square root of a rational function and gil/z) = g~^{z). The 
asymptotics of such determinant is computed in terms of multi-dimensional theta- 
functions associated to a hyperelliptic curve C of genus g > 1, which enter into 
the solution of a Riemann-Hilbert problem. Phase transitions for these systems 
are characterized by the branch points of C approaching the unit circle. In these 
circumstances the entropy diverges logarithmically. We also recover, as particular 
cases, the formulae for the entropy discovered by Jin and Korepin [13] for the XX 
model and Its, Jin and Korepin [12[ [TB] for the XY model. 
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1 Introduction 

One dimensional quantum spin chains were introduced by Lieb et. al. [I?] in 1961 as a 
model to study the magnetic properties of solids. Usually such systems depend on some 
parameter, e.g. the magnetic field. One of their most important features is that at zero 
temperature, when the system is in the ground state, as the number of spins tend to infinity 
they undergo a phase transition for a critical value of the parameter. As a consequence, 
the rate of the decay of correlation lengths changes suddenly from exponential to algebraic 
at the critical point. Furthermore, many examples of such chains are exactly solvable. 
Because of these reasons over the years the statistical mechanical properties of quantum 
spin chains have been investigated in great detail. 

More recently, Osterloh et al. [20], and Osborne and Nielsen [22] realized that the 
existence of non-local physical correlations at a phase transition is a manifestation of the 
entanglement among the constituent parts of the chain. Entangled quantum states are 
characterized by non-local correlations that cannot be described by classical mechanics. 
Such correlations play an important role in the transmission of quantum information. It 
is therefore essential to be able to quantify entanglement. In its full generality this is still 
an open problem. However, when a physical system is in a pure state and is bipartite, 
i.e. is made of two separate parts, say A and B, a suitable measure of the entanglement 
shared between the two constituents is the von Neumann entropy of either part [2] . In this 
situation the Hilbert space of the whole system is T^ab = T^a ^'H-q, where T^a and Tie are 
the Hilbert spaces associated to A and B respectively. Now, if pab is the density matrix 
of the composite system, then the reduced density matrices of A and B are 



Pa = tre Pab and = tiA Pab, 



(1.1) 
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where tiA and tie are partial traces over the degrees of freedom A and B respectively. The 
entropy of the entanglement of formation is 



In this paper we compute the entropy of entanglement of the ground state of a vast class 
of spin chains whose interaction among the constituent spins is non-local and translation 
invariant. These systems can be mapped into quadratic chains of fermionic operators by 
a suitable transformation and are generalizations of the XY model. We study the ground 
state of such systems, divide the chain in two halves and compute the von Neumann entropy 
in the thermodynamic limit of one of the two parts. If the ground state is not degenerate, 
then pab = |^g)(^g|- At the core of our derivation of the entropy of entanglement is the 
computation of determinants of Toeplitz matrices for a wide class of 2 x 2 matrix symbols. 
The explicit expressions for such determinants were not available in the literature. The 
appearance of Toeplitz matrices and their invariants in the study of lattice models is a 
simple consequence of the translation invariance of the interaction among the spins. Thus, 
Toeplitz determinants appear in the computations of many other physical quantities like 
spin-spin correlations or the probability of the emptiness of formation, not only the entropy 
of entanglement. Therefore, our results have consequences that go beyond the application 
to the study of bipartite entanglement that we discuss. 

Vidal et. al. [22] were the first to investigate the entanglement of formation of the 
ground state of spin chains by dividing them in two parts. The models they considered 
were the XX, XY and XXZ model. They computed numerically the von Neumann entropy 
of one half of the chain and discovered that at a phase transition it grows logarithmically 
with its length L. Jin and Korepin [Tl] computed the von Neumann entropy of the ground 
state of the XX model using the Fisher-Hartwig formula for Toeplitz determinants. They 
showed that at the phase transition the entropy grows like | logL, which is in agreement 
with the numerical observations of Vidal et. al. For lattice systems that have a conformal 
field theory associated to it the logarithmic growth of the entropy was first discovered by 
Holzhey et. al. [10] in 1994. This approach was later developed by Korepin [15], and by 
Calabrese and Cardy [1]. Its, Jin and Korepin [T2l [T3] determined the entropy for the 
XY model by computing an explicit formula for the asymptotics of the determinant of a 
block- Toeplitz matrix. They expressed the entropy of entanglement in terms of an integral 
of Jacobi theta functions. 

Consider a. p x p matrix- valued function on the unit circle S: 



trpAlogpA = S{pb) 



trpB logpB 



(1.2) 



oo 



z\ = 1. 



fc=— oo 



A block- Toeplitz matrix with symbol is defined by 



= {fj~k)o<j,k<L-l- 
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Furthermore, we shall denote its determinant hj Dl = detT^ The main ingredient 
of the computation of Its, Jin and Korepin was to use the Riemann-Hilbert approach to 
derive an asymptotic formula for the Fredholm determinant 

Dl(A) = detrz.M =det(J-KL), (1.3) 

where K^, is an appropriate integral operator on L^(H,C^). The symbol of the Toeplitz 
matrix Tl[lp] was 



where 

acos6' — 1 — Z7asin0 

9{0) 



I a cos 9 — 1 — i'ja sin 9 \ 

Keating and Mezzadri [T8l [T9] introduced families of spin chains that are characterized 
by the symmetries of the spin-spin interaction. The entropy of entanglement of the ground 
state of these systems, as well as other thermodynamical quantities like the spin-spin 
correlation function, can be determined by computing averages over the classical compact 
groups, which in turn means computing determinants of Toeplitz matrices or of sums of 
Hankel matrices. These models are solvable and can be mapped into a quadratic chain of 
Fermi operators via the Jordan- Wigner transformations. One of the main features of these 
families is that symmetries of the interaction can be put in one to one correspondence with 
the structure of the invariant measure of the group to be averaged over. If the Hamiltonian 
is translation invariant and the interaction is isotropic, then the relevant group over is U(A^) 
equipped with Haar measure. In turn such averages are equivalent to Toeplitz determinants 
with a scalar symbol. These systems are generalizations of the XX model. 

In this paper we consider spin chains whose interaction is translation invariant but the 
Hamiltonian is not isotropic. These are generalization of the XY model. The Fredholm 
determinant that we need to compute has the same structure as (11.31) . but now the 2x2 
matrix symbol is 

*W-( "i^'V (1-5) 



where function g{z) is defined by 



g{z) := J , (1.6) 



and "pi^z) is a polynomial of degree In. We recover the XY model if we set 

a(l -7) 2 , «(1 +7) (. 7^ 

p{z) = z -z-\ . (1.7) 



In the above equation a = 2/h, where h magnetic field, and 7 measures the anisotropy of 
the Hamiltonian in the XY plane. 
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Rc(A) 



Figure 1: The contour r(e) of the integral in equation (12.1 1) . The bold lines (— oo, —1 — e) 
and (1 + e, oo) are the cuts of the integrand e(l + e, A). The zeros of Dl{\) are located on 
the bold line (—1, 1). 

2 Statement of results 

Following [13] and [13] , we will identify the limiting von Neumann entropy for the systems 
that we study with the double limit 



In the above formula r(e) is the contour in figure [T], -Dl(A) is the determinant of the 
block- Toeplitz matrix T;,[<I>] with symbol fll.Sp and 



The explicit Hamiltonians for the family of spin systems that we consider and their con- 
nection to formula (12 .ip will be discussed in detail in sections [3] and HI 

One of the main objectives of this paper is to compute the double limit (12. ip . which, as 
we shall see, can be expressed as an integral of multi-dimensional theta functions defined 
on Riemann surfaces. Thus, in order to state our main results, we need to introduce some 
definitions and notation. 

Let us rewrite the function (ll.6p as 



S{pa) = lim lim i — 




(f e(l + e,A);^log(Di(A)(A2-l)-^)dAj . (2.1) 




(2.2) 



2n 



= n 



(2.3) 
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Re(A,) 


< 


Re(Aj), i<j 






Im(A,) 


< 


Im(Aj), i<j, 1 


Ail, < 1, Re(Ai) 


= Re{\j 


Im(A,) 


< 


Im(Aj), i>j, 1 


Ai , \Xj\ > 1, Re(Ai) 


= Re{\j 



where the zj^s are the 2n roots of the polynomial p{z). This representation of g{z) will be 
used throughout the paper. We fix the branch of g{z) by requiring that g{oo) > on the 
real axis. The function g{z) have jump discontinuities on the complex 2;-plane. In order to 
define its branch cuts we need to introduce an ordering of the roots Zj. Let 

} (2.4) 

where the above is merely an equality between sets, and we do not necessarily have, for 
example, Aj = Zi. We order the Aj's such that 



(2.5) 



This ordering need not coincide with the ordering Zj-s,. If necessary, we can always assume 
that one of the 2;"^ has the smallest real part and set Ai = z~^ . This choice is equivalent 
to taking the transpose of Ti[$]. The branch cuts for (7(2) are defined by the intervals Sj 
joining A2i-i and X^o 

Si = [Asi-i, A2i], i = l,...,2n. (2.6) 
Therefore, gi^z^ has the following jump discontinuities: 

g^{z) = -g.{z), ^gS„ (2.7) 

where g±i^z^ are the boundary values of (7(2) on the left/right hand side of the branch cut. 
Now, let £ be the hyperelliptic curve 

L:w' = ^{z-X^. (2.8) 

i=\ 

The genus of £ is = 2n — 1. We now choose a canonical basis for the cycles {oj, 6j} on £ 
as shown in figure [2], and define dcu^ to be 1-forms dual to this basis, i.e. 

dwj = / dwj = Iljj. (2.9) 

Furthermore, let us define the g ^ g matrix 11 by setting (n)jj = Iljj. The theta function 
d : — > C associated to L is defined by 

e{-t) := J2 e^-"'-n"'+2..^-^. (2.10) 



while the theta function with characteristics ~t and 6 is defined by 

e 



e 

t 



exp{2^n(^^^^ + ^-t■^ + \-t■^)U^r + ^ + U^ 1(2.11) 



where ~t and 6 are (yf-dimensional complex vectors. 

Our main results are summarised by the following two theorems. 
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\z\ = 1 



A] 




Figure 2: The choice of cycles on the hypereUiptic curve C. The arrows denote the orien- 
tations of the cycles and branch cuts. Note that we have Ai = ^f^. 



Theorem 1. Let be the Hamiltonian of the one- dimensional quantum spin system 
defined in equation ^3.11\) . Let A he the subsystem made of the first L spins and B the one 
formed by the remaining M — L. We also assume that the system is in a non- degenerate 
ground state |\E'g) and that the thermodynamic limit, i.e. M — > oo, has been already taken. 
Then, the limiting (as L oo) von Neumann entropy h2. 1\) is 



where ~^ is a 2n — 1 vector whose last n entries are 1 and the first n — 1 entries are 0. 

The parameter r in the argument of 9 is introduced in section [6] and is defined in 
equation fl6.11l) . while the expression of [3{\) is 



Theorem [T] generalizes the result by Its et al. [T2l [T3] for the XY model. In that case 
the genus of the of C is one, and the theta function in the integral reduces to the Jacobi 
theta function 9^. However, for the XY model the integral fl2.12p can be expressed in term 
of the infinite series 




(2.12) 



A + 1 

A- r 



(2.13) 




(2.14) 



where the numbers /i^ are the solutions of the equation 




(2.15) 
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and 0" is or 1 depending on the strength of the magnetic field. The zeros of the one 
dimensional theta function are all known, so that the numbers /^m can be described by the 
explicit formula 



Moreover, as it was shown by Peschel [2T] (who also suggested an alternative heuristic 
derivation of equation (12.141) based on the work of Calabrese and Cardy [3]), the series 
(12.141) can be summed up to an elementary function of the complete elliptic integrals 
corresponding to the modular parameter r. 

It is an open problem whether an analogous representation of the integral (12.121) exists 
for 5' > 1. 

The next step consists of understanding what happens to formula (I2.12p when we 
approach a phase transition. The hyperelliptic curve C, and hence all the parameters in 
the integral (I2.12p . are determined by the roots of the polynomial p{z) which defines the 
symbol (11.51) . In section [3] we discuss how the coefficients of p(z) are related to the the 
Hamiltonians of the spin chains. In the case of the XY model p{z) is given by equation (11.71) : 
since the degree of p{z) is two the roots can be easily determined as a function of 
the parameters a and 7. It was shown by Calabrese and Cardy [1] that when a = 1 
— or the magnetic field h = 2 — the XY model undergoes a phase transition and the 
entropy diverges. Jin and Korepin [13] showed that when 7 approaches 0, i.e. the XY 
model approaches the XX model, and a < 1, then the entanglement entropy diverges 
logarithmically. Its et. al. [121 US] discovered that the divergence of the entropy for the 
XY and XX model corresponds to the roots (12. 4p of (12.81) approaching the unit circle. 

This phenomenon extends to the family of systems that we study. In other words, 
a phase transition manifests itself when pairs of roots of (12. 8 p approach the unit circle; 
one root in each pair is inside the unit circle, the other outside. As we shall see, in these 
circumstances the entropy of entanglement diverges logarithmically. From (12. 4p we see that 
if Xj is a root of (12. 8p so is Xj^. Moreover, since (12. 8p is a polynomial with real coefficients, 

if Xj is complex then Xj and A^- ^ will be roots of (12. 8p too (see figure [3]). Now, suppose that 

Xj approaches the unit circle and \Xj\ < 1, then | Aj| ^ > 1 and A^- ^ will also be approaching 
the unit circle with 



At a phase transition the behavior of the entropy of entanglement is captured by 





0. 



Theorem 2. Let the m pairs of roots Xj, Xj , j = 1, . . . ,m, approach together towards 

the unit circle such that the limiting values of Xj, Xj ^ are distinct from those of X^, A^/ if 
j k, then the entanglement entropy is asymptotic to 




m 



j = 1, . . . ,m. 



(2.16) 
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A. / 



Figure 3: The location of one of the roots (12.41) . say \j determines the positions of other 
three: Aj, 1/Aj and 1/Aj 

From the integral (12.1 1) it is evident that in order to prove theorems [T] and [2] we need an 
explicit asymptotic formula for the determinant Dl{\). Indeed, the following proposition 
gives us an asymptotic representation for the determinants of block- Toeplitz matrices whose 
symbols belong to the family defined in equations fll.Sp and (11.61) . 

Proposition 1. Let be the set 

n,:={XeR: |A| > 1 + e}. (2.17) 

Then the Toeplitz determinant D admits the following asymptotic representation, which 
is uniform in \ E 

, 9xf^(/5(A)^ + i)^(/5(A)e'-i) / , j.\ , ^ 

D^{\) = {\-\^f^^ ^(l + 0(p-^)), L-.00, (2.18) 

Eere p is any real number satisfying the inequality 

1 < p < min{|Aj| : |Aj| > 1}. 

Remark 1. The first factor in the right hand side of equation 1^2. 18\) corresponds to the 
"trivial" factor, G[^] of the general Widom's formula liS. which we discuss in detail in 
section [31 while the ratio of the theta functions provides an explicit expression of the most 
interesting part of the formula — Widom's pre-factor E[^] = det (Too[$]Too[$~^]), which 
is given in formula l^5.2\) . 
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Remark 2. The Asymptotic representation 1^2. 18\) is actually valid in a much wider domain 
of the complex plane A. Indeed, it is true everywhere away from the zeros of the right hand 
side, which, unfortunately, in the case of the genus g > 1 is very difficult to express in a 
simple closed form — one faces a very transcendental object, i.e. the theta- divisor. This 
constitutes an important difference between the general case and that one with g = 1 studied 
in fW^ and flS^ . where the zeros of equation Ii2.15\) can be easily evaluated. 



3 Quantum spin chains with anisotropic Hamiltoni- 
ans 

The XY model is a spin-1/2 ferromagnetic chain with an exchange couphng a in a constant 
transversal magnetic field h. The Hamiltonian is H = hHa with Ha given by 

M-1 M-1 

= -f E + ^)^J^]Vi + (1 - 7)a|a|^i] - E (3.1) 

j=0 j=0 

where {a^,a^,a^} are the Pauli matrices. The parameter 7 lies in the interval [0, 1] and 
measures the anisotropy of Ha- When 7 = (13.11) becomes the Hamiltonian of the XX 
model. In the limit M —>■ 00 the XY model undergoes a phase transition at = 1. 

It is well known that the Hamiltonian (13.11) can be mapped into a quadratic form of 
Fermi operators and then diagonalized. To this purpose, we introduce the Jordan- Wigner 
transformations. Let us define 



m2i+i 



The inverse relations are 



HaMaf and m,i = ma]) af . (3.2) 
d=o J \j=o / 



'i-i 



= I JJ^"^2j"i2i+l ) "^2Z+l, 

i=o / 
Z-1 \ 

JJim2jm2j+i m2i (3.3) 



These operators obey the commutation relations {mj,mk} = 25 jk but are not quite Fermi 
operator since they are Hermitian. Thus, we define 

hi = (m2/+i - im2i)/2 and b] = {m2i+i + im2i)/2, 

which are proper Fermi operator as 

{bj, bk} = and {bj, fej,} = Sjk- 
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In terms of the operators 6j's the Hamiltonian (13.11) become^ 

Ha 



a 
2 



A/-1 



j=0 



M-1 

i=o 



(3.4) 



It turns out that the expectation values of the operators (13.21) with respect to the ground 
state |\l/g) are 



0, 



^jk + 'i(yCM)jk, 



(3.5) 
(3.6) 



where the correlation matrix Cm has the block structure 



C 



M 



( Cu 

C21 



Ml 



with 



jk 



C12 
Cm2 



-9k^j 



ClM \ 
C2M 

Cmm / 



(3.7) 



9j~k 




For large M, the real numbers gi are the Fourier coefficients of 

acos^ — 1 — 27asin6' 
= I Q — \ '■ 

In other words, Cm is a block- Toeplitz matrix with symbol 



g{e) 

-g-\9) 



(3.^ 



(We outline the derivations of formulae (13.51) and (13.61) for the family of systems (13.101) that 
we study in the appendices B and C.) 

Equation (13.51) is a straightforward consequence of the invariance of under the map 
bj (— s> —bj] for the same reason the expectation value of the product of an odd number of 



mj's must be zero. Formula (13. 6p was derived for the first time by Lieb et al. [T7]. The 
expectation values of the product of an even number of the mj's can be computed using 
Wick's theorem: 



(^'„| rrii.m 



m 



'j2n l^g) = E^ {^^Y W (contraction of the pair) 

all pairings all pairs 



(3.9) 



^This is strictly true only for open-end Hamiltonians. If we impose periodic boundary conditions, then 



the term b\^_^ba in (|3.4p should be replaced by J| 



2h\h 



^l/-i^o- However, because we are 



interested in the limit Af ^ 00, the extra factor in front of b\j_^bQ can be neglected. 
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where a contraction of a pair is defined by (^'g| mj^mj^ ^i-nd p is the signature of the 
permutation, for a given pairing, necessary to bring operators of the same pair next to 
one other from the original order. Many important physical quantities, including the von 
Neuamnn entropy and the spin-spin correlation functions, are expressed in terms of the 
expectation values fl3.9p . 

In this paper we study generalizations of the Hamiltonian (13.41) that are quadratic in 
the Fermi operators and translation invariant. More explicitly, we consider the family of 
systems 

M-l 

-2Y,h]h, (3.10) 

i=o 

with cyclic boundary conditions. In terms of Pauli operators this Hamiltonian becomes 



Ha 



a 



'M-l 



a \^ 



jk 



0<j<k<AI-l 



y 

J ^k 



k-l 

n 



cr, 




(3.11) 



The translation invariance of the interaction implies that Ajk = Aj^k and Bjk = Bj^k, 
and the cyclic boundary conditions force A and B to be circulant matrices. Furthermore, 
since Ha is a Hermitian operator, the matrices A and B must be symmetric and anti- 
symmetric respectively. Now, let us introduce two real functions, 

a : Z/MZ — >R and b : Z/MZ — > M, 

such that 

a{j -k) = aAj_k - 26 jk and b{j - k) = aBj_k, j, k e Z/MZ. (3.12) 

Since A is symmetric and B anti-symmetric, we must have 

= «(j) and b{-j) = -b{j). 

We shall consider systems with finite range interaction, which implies that there exists a 
fixed n < M such that 

a{j)=b{j) = for j>n. (3.13) 

In the appendices B and C we derive the expectation values in the ground state of 
the Jordan- Wigner operators m^' s. They have the same structure as the expectation 
values (13.51) and (13.61) . but now in the limit as M — * cxd the symbol (13. 8p of the correlation 
matrix Cm is replaced by 











(3.14) 
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where 

n 

qiz) = {a{j) - ih{j)) (3.16) 

j=-n 

p{z) = z^q{z). (3.17) 

4 The von Neumann entropy and block- Toeplitz de- 
terminants 

We now concentrate our attention to study the entanglement of formation of the ground 
state l^'g) of the family of Hamiltonians (13.101) . Since the ground state is not degenerate, 
the density matrix is simply the projection operator |^'g)(^'g|. We then divide the system 
into two subchains: the first one A containing L spins; the second one B, made of the 
remaining M — L. We shall further assume that 1 <^ L <C M . This division creates a 
bipartite system. The Hilbert space of the whole system is the direct product TYab = 
7i A ® 'Hb , where Ha and Hb are spanned by the vectors 

n(&])'M*vac) and l[{b]Y^\^..c), r, =0,1, 

respectively. The vector l^'vac) is the vacuum state, which is defined by 

&il*vac)=0, J=0,...,M-1. 

Our goal is to determine the asymptotic behavior for large L, with L = o(M), of the von 
Neumann entropy 

5'(pa) = -trpAlogpA, (4.1) 

where pa = trePAB and pab = l*g)(*gl- 

It turns out that after computing the partial trace of pab over the degrees of freedom of 
B, the reduced density matrix pA can be expressed in terms of first L Fermi operators that 
generate a basis spanning Ha- As a consequence, only the submatrix Cl formed by the first 
2L rows and columns of the correlation matrix (13.71) will be relevant in the computation of 
the entropy (14.11) . Now, Cl is even dimensional and skew-symmetric. Furthermore, since 



g (e-^) = g (e^^) 

its Fourier coefficients are real, therefore there exists an orthogonal matrix V that block- 
diagonalizes Cl. 



VC,V = ^.,(°^ I), (4,2) 

j=0 ^ / 
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where the ±zz/jS are imaginary numbers and are the eigenvalues of the block- Toephtz matrix 
Cl = TlI'^], where ip is the symbol (13.141) . 
Let us introduce the operators 

Cj = (dsi+i - id2j)/2, j = 0, . . . , L - 1, (4.3) 

where 

2L-1 

dj = ^ikrrik. (4.4) 

k=Q 

Since V is orthogonal {dj^dk} = 26j and the cjs are Fermi operators. Combining equa- 
tions (14.21) . (14.31) and (14. 4p . we obtain the expectation values 

(*g|c,|*g) = (*g| c,-Cfc I *g) =0, (4.5) 

(*g|c]cfc|*g) = 6,k^-^. (4.6) 

The reduced density matrix pA can be computed directly from these expectation values. 
We report this computation in appendix A. We have 



J=0 



In other words, as equations (14.51) and (14.61) already suggest, these fermionic modes are in 
a product of uncorrelated states, therefore the density matrix is the direct product 

PA = (g)Pj With p^ = —^c]c, + ^^c,c]. (4.8) 

j=0 

Since (1 + i'j)/2 and (1 — ^'j)/2 are eigenvalues of density matrices they must lie in the 
interval (0, 1), therefore, 

-l<iy^<l, j = 0,...,L-l. 

At this point the entropy of the entanglement between the two subsystems can be easily 
derived from equation (14. ip : 

L-l 

^(PA) = 5^e(l,i.,), (4.9) 

j=0 

where e{x, u) is defined in equation (12. 2p . Using the residue theorem, formula (14.91) can be 
rewritten as 

SM ^ ln.^£^^((-l)^|3^)e(l + ,^ 

= hm^/ eil + e,X)^^^^^^dX (4.10) 
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where r(e) is the contour in figure [Hand 

L-l 

^L(A)=(-i)"n(^'-^i) (4.11) 

is the determinant of the block- Toephtz matrix Tl[$](A) with symbol (11. 5p . 

The integral (I4.10p was introduced for the first time by Jin and Korepin [T3] to com- 
pute the entropy of entanglement in the XX model. In this case g~^{9) = g{6) and Dl{X) 
becomes the determinant of a Toeplitz matrix with a scalar symbol. Keating and Mez- 
zadri [181 CH] generalized it to lattice models where Dl{\) becomes an average over one 
of the classical compact groups. Its et al. [Ill [13] computed the same integral for the XY 
model, for which Dl{X) is the determinant of a block- Toeplitz matrix with symbol (11.41) . 
Following the same approach of Its et al., in this paper we express Dl{X) as a Fredholm 
determinant of an integrable operator on Iv^(S, C^) and solve the Riemann-Hilbert problem 
associated to it. This will give an explicit formula for Dl{X), which can then be used to 
compute the integral (14.101) . 



5 The Asymptotics of Block Toeplitz Determinants. 
Widom's Theorem 

A generalization of the strong Szego's theorem to determinants of block- Toeplitz matrices 
was first discovered by Widom p^l25]. Consider a.pxp matrix symbol ip and assume that 

oo / oo \ 1/2 

\\'^\\= llv^fcll+l 5Z I^NIv^fcllM < oo- 

fc=— oo \fe=— oo / 

The norm that appear in the right-hand side of this equation is the Hilbert-Schmidt norm 
of the p X p matrices that occur. In addition, we shall require that 

detip{z)y^O and A||2|=i argdet ^9(2) = 0. 

Widom showed that if one defines 

:=exp(^^^logdet(^(2;)^^ (5.1) 

then 

E[cp] := ^im = det {T^[cp]T^[^-']) , (5.2) 

where Too[v5] is a semi-infinite Toeplitz matrix acting on the Hilbert space of semi-infinite 
sequence of ]?- vectors: 



k=0 



VfcGC^ ^||vfc||2 < 00 I . 



A;=0 
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Formulae (15.11) and (15.21) reduce to Szego's strong limit theorem when p = 1. Although 
this beautiful formula is very general, it is difficult to extract information from the right- 
hand side of equation (15.21) and determine formulae that can be used in the applications. 
The advantage of our approach is precisely to derive explicit formula for the leading order 
term of the asymptotics of block- Toeplitz determinants whose symbols ^{z) belong to the 
one-parameter family defined in (11.51) . 

A starting point of our analysis is the asymptotic representation of the logarithmic 
derivative (with respect to the parameter A) of the determinant -Dl(A) = detTL[$](A) 
in terms of 2 x 2 matrix-valued functions, denoted by U±{z) and V±{z), which solve the 
following Wiener-Hopf factorization problem: 

^z) = U+{z)U.{z) = V.{z)V+iz), (5.3) 

U^{z) and V_{z) {U^{z) and V+{z)) are analytic outside (inside) 

the unit circle S, (5.4) 

f/_(cx)) = 1/_(oo)=/. (5.5) 

Now, let us fix e > and define the set 

a := {A G M : |A| > 1 + e}. (5.6) 

In the next section we will show that for every A G fi^ the solution of the above 
Wiener-Hopf factorization problem exists, and the corresponding matrix functions, U±{z) 
and V±{z) satisfy the following uniform estimate: 

\u-{z)\, \v-{z)\ < a, yzev±, vAGa, (5.7) 

where the notation P+ (^^-) is used for the interior (exterior) of the unit circle S. Moreover, 
generalizing the approach of [121 [IS] we will obtain the multidimensional theta function 
explicit formulae for the functions U±{z) and V±{z). 

The asymptotic representation of the logarithmic derivative dlogDL(A)/dA is given by 
the following theorem: 

Theorem 3. Let A G ^2^, and fix a positive number R > 0. Then, we have the follow- 
ing asymptotic representation for the logarithmic derivative of the determinant -Dl(A) = 
detTi[$].- 

^log/^.(A) ^ 



2" 



dX ^ ' ' 1-A 



dz 



{U'^{z)U~\z) + V;\z)V[{z))^-\z) 

+rL(A), (5.8) 
where (') means the derivative with respect to z, the error term tl^X) satisfies the estimate 
\rL{\)\<Cp-^, A G a n {|A| < i?}, L>1, (5.9) 
and p is any real number such that 1 < p < min{|Aj| : \Xj\ > 1}. 



18 



A. R. Its, F. Mezzadh and M. Y. Mo 



This theorem, without the error term estimate is a specification of one of the classical 
results of H. Widom [21] for the case of the matrix generators ^{z) whose dependence on 
the extra parameter A is given by the equation 

^{z) =^{z-X) =zA/ + $(^;0). 

The estimate (15.91) of the error term as well as an alternative proof of the theorem itself 
in the case of curves of genus one is given in [12] and [13]. The method of [12] and [13] 
is based on the Riemann-Hilbert approach to the Toeplitz determinants [3] and on the 
theory of the integrable Fredholm operators [HI [9]; its extension to symbols (11. 5p . where 
the polynomial p{z) entering in (II. 6p is of arbitrary degree, is straightforward. Indeed, the 
following generalization of theorem [3] follows directly from the analytic considerations of 

m. 



Theorem 4. Suppose that the matrix generator $(2;) is analytic in the annulus, 

Vs = {l-5 <\z\< 1 + 5}. 

Suppose also that $(2) depends analytically on an extra parameter fi and that it admits a 
Wiener- Hopf factorisation for all fi from a certain set Ai. Finally, we shall assume that 
the matrix functions 

^1^, U^{z), andV^iz) 

are uniformly bounded for all ^ & Ai and all z from the respective domains, i.e. Vs in the 
case of^{z), ^~^{z), and d^{z)/dfi, and V± in the case of U±{z) andV±{z). Then, the 
logarithmic derivative of the determinant Dl{ij) = det Tl[^] has the following asymptotic 
representation: 

^logDM = _ / tr $-1(2;)— ^ — + — / tr (<l>-iy(z) — 



dfi 27ii Jh V dfi J z 2711 \ dfi 

+^^tr (u:,{z)U^\z)^<l>-\z) + V-\z)Vi{z)<l>-\z)^^ dz 
+rM, (5.10) 
where the error term ri(/i) satisfies the uniform estimate 

|rL(/i)| < Cp-^, fieM, L>1, (5.11) 
and p is any positive number such that 1 < p < 1 + 5. 



This theorem, without the estimate of the error term and with much weaker assumptions 
on the generator ^{z), is exactly the classical result of Widom from |24j . 
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Remark 3. Denote 

u±{z) = V^-^i^z), and v±{z) = U^-^^z) 

so that 

^~^{z) = u+{z)u^{z) = v^{z)v+{z). 
Then, equation 1^5. 10\) can be re-written in a more compact way: 

d^{z)\ dz 



+rM- (5.12) 
This form in which this result is formulated in l2^ . 

Theorem H] can be used to strengthen the statement of theorem [3] by removing the 
dependence of the constant C on i? in the estimate (15.91) . This leads to the following 
extension of theorem [S] 

Theorem 5. Let VL^ he the set defined in \5. 0\] and let $(z) he the symhol defined in i\1.5\) . 
Then we have the following asymptotic representation of the logarithmic derivative of the 
determinant Dl{X) = det Tl[^] for all A G fi^.' 

^logD^(A) = -^^L + l-J^tT[{U'4z)U~\z) + V-\z)Vi{z))<!>-\z) 

+rL(A), (5.13) 

where (') means the derivative with respect to z, the error term tl^X) satisfies the uniform 
estimate 

|r•L(A)|<^p-^ AG a, L>1 (5.14) 
and p is any real numher such that 1 < p < min{|Aj| : |Aj| > 1}. 

Proof. Let i? > 1 + e and denote Ci the constant C from estimate fl5.9p . Take now A G ^2^, 
|A| > R and set 

1 . . f ^ , , 1 1 
p = -eM = <peR: \p\< — < 



dz 



A I '"'-i? 1 

By trivial algebra, we arrive at 

detZ}z.(A) = (-A2)^det^i(/i), 

where D^^p) = detTi[<l'] and 

1 



ipg[z) 
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From this relation it also follows that 

^\ogdet Dl{\) = ^ - ^^logdetDi(/i), (5.16) 

and hence the asymptotic analysis of the logarithmic derivative dlogdet Z)i;,(A)/dA for 
|A| > -R is reduced to that one of the logarithmic derivative dlogdet Di(/i)/d/i for ji G 

M^{^^e^■.\^^\<^<^^]. 

Firstly, we notice that for all yU G and z E Vs the functions ^{z), and 
d^{z)/dn are uniformly bounded. Secondly, we have that 

^z) = ^^z) = \u4z)U4z) = ^V4z)V4z), 

and hence the matrix valued functions U±{z) and V±{z) defined by the relations 

[7+(z) = ^f/+(z), V+{z) = ^V+{z), U.{z) = U.{z), V.{z) = V.{z) 

provide the Wiener- Hopf factorization of the generator $(2). Moreover, because of the 
estimates fl5.7l) . the functions U±{z) and V±{z) are uniformly bounded for aX\ n E M. and 
z G Hence, all the conditions of theorem H] are met, and we can claim the uniform 
asymptotic representation (15.101) of the logarithmic derivative of the determinant Dl[^) 
with the symbols $, [/, and V replaced by $, U and V respectively. We shall also use the 
notation tl^^) and C2 for the error term and constant C from the corresponding estimate 
(15.111) respectively. 

The specific form (I5.15P of dependence of the generator ^{z) on the parameter /i implies 
that 

i-w^ = _L_ ( . -sW), (5.17) 



and 



Hence 



.i,,,,9$(2)_ 1 -^,g~\z)g'{z) 



(<l>-^)'(^)— ^ = );^' (5.18) 



, 1, ^9<l>(z)\ 2/i 2A 



dn ) 1 — /i^ 1 — A 



2 



and equation (15.101) for the determinant -Dl(/^) becomes 
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+rM, (5.19) 

with 

\fM\<C2P^'', /ieTW, L>1. (5.20) 
Observe now that equation fl5.17p can be rewritten as 

This relation, together with the obvious fact that 

Kiz)U-\z) = U'4z)U~\z) and V-\z)Viiz) = V~\z)Viiz), 
allows to transform fl5.19p into the asymptotic formula 

— log Dl{Ij,) 



^ tr [{U'4z)U-\z) + V-\z)Vi{z)) ^-\z)] dz 
+rM. (5.21) 

The substitution of this relation into the right hand side of equation f l5.16p yields the 
following asymptotic formula — which is complementary to the equation (15.81) — 



dz 



{U'^{z)U-\z) + V-\z)Vi{z))<l>-\z)_ 
+rL(A), (5.22) 
with the error term rL{X) satisfying the estimate 

\rLW\<j^p-\ Aean{|A| >/?}, L>1. (5.23) 

Choosing 

C = max{Cii?,C2}, 

we arrive at the statement of the theorem, but with a better estimate for the error term 
rL(A) than that one in fl5.14p . 

In order to improve the estimate fl5.23p . we notice that since ^{z) becomes the identity 
matrix as /i — *• 0, the Wiener-Hopf factorization of ^{z) exists for all fi from the small 
complex neighbourhood 

Mo = {fieC:\ii\<eo< -} 



22 



A. R. Its, F. Mezzadh and M. Y. Mo 



of the point /i = 0. In particular, this imphes that the Wiener- Hopf factors, U±{z) and 
V±{z), admit an analytic continuation to the disc A^o and that the validity of the formulae 
(15.191) and fl5.20p can be extended to the set 

Moreover, from equation fl5.19p it follows that r2,(/i) is analytic in the disc M.q and that 
^1,(0) = 0. In order to see that the latter equality is true, one has to take into account that 
U±{z) = V±{z) = I for all z and fi = and the evenness of Di{fi) as a function of fi. Now, 
define 

rM = . 

The function "rilfi) is analytic in the disc JAq and satisfies the estimate (15.201) uniformly 
for n & C^i = {|/i| = e'} and for any < e' < eo. With the help of the Cauchy formula, 

we conclude that 

\fM\<Cp^, lAi|<eo/3, L>1 

or 

|fi(/i)| < C|/i|p^, |/i|<eo/3, L>1. 
The last inequality combined with (15.201) allows to replace it by the estimate 

< C\fi\p^, fieM, L>1, 

which, in turn, transforms estimate (15.231) into the estimate 



rL(A)|<^p"^, A G an{|A| > i?}, L>1, (5.24) 



lA 



and hence yields the correction term as announced in (15.141) . This completes the proof of 
the theorem. □ 



6 The Wiener-Hopf factorization of ^{z) 

In this section we will compute the Wiener-Hopf factorization of ^{z). We will express the 
solution in terms of theta functions on a hyperelliptic curve C 
From the equality 



{l-X^)a3^-\z)a3 = ^z), (T3 = 
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we can express V in terms of U as follows: 

V.{z) = a^Uz'a, 

V+{z) = asU-\z)asil-\''), A ^ ±1. (6.1) 

Therefore, we only need to compute U{z). To do so, first note that $(2;) can be diagonalized 
by the matrix 

Q{z)=('^'^ -'^'A. (6.2) 



Indeed, it is straightforward to see that 



4(2) = Q(j)AQ-'(:) 
^ .,A+1 



A- 1 

The function Q{z) has the following jump discontinuities on the z-plane: 



Q+{z) = Q-{z)ai, 2 e Si 

1 

1 



(71 



where the branch cuts Sj are defined in (12. 4p . (12. 5p and (12.61) and Q±{z) are the boundary 
values of Q{z) to the left /right of Sj. It also has square-root singularities at each branch 
point with the following behavior: 



Qiz) = Q±^iz) - Z^Z 



± 
i 1 



where Q±i{z) are functions that are holomorphic and invertible at zf. 
Let us define 

S{z) = U4z)Q{z)A-\ \z\ > 1, 

S{z) = U+{z)-'Qiz), \z\ < 1. (6.3) 

By direct computation we see S{z) is the unique solution of the following Riemann-Hilbert 
problem: 

S+{z) = S-{z)ai, 2 G Sj, i = l,...,n 

S+{z) = S4z)AaiA~\ zeE„ i = n+l,...,2n (6.4) 
lim S{z) = Q{oo)A-\ 

where, as before, S±{z) denotes the boundary values of S{z) to the left/right of the branch 
cuts. The matrix function S{z) is holomorphic and invertible everywhere, except on the 
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cuts Sj, where it has the jump discontinuities given in (16.41) . and in proximity of the branch 
points, where it behaves hke 



[z - z^^)^'2 \ / 1 -1 



S{z) =S^,{z)^^- 2> ^j, z^zf, |z,|<l, (6.5) 

S{z) = S^,{z) (^^' - 'p^^' -^)a-\ z-.zf, \z^>l. 

where S±i{z) are holomorphic and invertible at zf. 

The Riemann-Hilbert problem (16.41) can be solved in terms of the multi-dimensional 
theta functions (I2.10p . However, before we compute explicitly S{z), we need to introduce 
further notions and properties of 9. 

Throughout the rest of this section we shall use the definitions (12.81) of the hyperelliptic 
curve C and (12.101) of the theta function associated to C Furthermore, recall that the 
choice of the canonical basis for the cycles is described in figure [2] and that the normalized 
1-forms dual to this basis are defined in equation (12.91) . Let us introduce some basic 
properties of the theta functions. The proofs of such properties can be found in many 
standard textbooks in Riemann surfaces like, for example, [H]. 

Proposition 2. The theta function is quasi-periodic with the following properties: 

9{-t + = ^(T), (6.6) 



9{^ + U'M) = exp 2ni (Pl,^^ - (li,^!!"^^ 



(6.7) 



where (■, ■) denotes the usual inner product in 

A divisor D of degree m on a hyperelliptic curve £ is a formal sum of m points on C, 

i.e. 



m 



=1 



Let us introduce the Abel map ui : C — > by setting 



CP rp 

I po Jpo 



where po is a chosen base point on C and Ui are the normalized 1-forms given in (12. 9p . In 
what follows we shall set po = = -^i- The composition of the theta function with the 
Abel map has g zeros on S. The following lemma tells us where the zeros are. 

Lemma 1. Let D = J2i=i be a divisor of degree g on C, then the multivalued function 



9{u{p) -uj{D) - K) 
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has precisely g zeros located at the points di, i = 1, . . . , g . The vector K = {Ki, . . . , Kg) is 
the Riemann constant 

The hyperelliptic curve C can be thought of as a branched cover of the Riemann sphere 
C U {oo}. Indeed, a point p E C can be identified by two complex variables, p = {z,w), 
where w and z are related by equation (12.81) . We shall denote by Ci the the Riemann 
sheet where g{oo) > on the real axis, and by C2 the other Riemann sheet in C Thus, a 
function f on L can be thought of as a function in two complex variables: 

f{p) = f{z,w). 

Consider the map 

T : C/ U^^i — >C 
T{z) = {z,w), 

where the branch of w is chosen such that {z, w) is on Ci. A function f on £ then defines 
the function / o T on C/ U^!!^ by 

foTiz) = f{z,w). 

For the sake of simplicity, and when there is no ambiguity, we shall write f{z) instead of 
/ o T{z) and f{p) instead of f{z,w). 

Abelian integrals on C can be represented as integrals on the Riemann sheet with jump 
discontinuities. To do so, let us first define a Jordan arc S as in figure HI Let f{z, w) be a 
function on jC and f{z) = / o T{z). Then an Abelian integral on C, 

i{p)= ['fipW: 

defines the following integral on C: 

I{z)= r foT{z')dz\ 
Jxi 

where the path of the integration does not intersect S/{Ai}. Such integral will in general 
have jump discontinuities along E, and its value on the left hand side of S will be denoted 
by I{z)+, while its value on the right hand side of S will be denoted by I{z)-. 
Let p be the hyperelliptic involution that interchanges the two sheets of £, i.e. 



p{z,w) = (2, -w). 
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The action of p on /(z) is given by 

p{f){z) = f{z,-w) (6.8) 

i.e it is the function evaluated on C2. Similarly, the action of p on an integral I{z) is 
defined by 

P(/)W= fpUWW (6.9) 




Figure 4: The Jordan arc S connects all the branch points and extends to infinity on the 
left hand side of Ai and on the right hand side of A4„. All branch cuts belong to S and are 
denoted by Sj, while the intervals between the branch cuts are denoted by Ej. 



From proposition [2] we see that the composition of the Abel map uj with 9 has the 
following jump discontinuities when considered as a function on C: 

Lemma 2. Let z he a point on C, and let H he a Jordan arc joining all the hranch cuts as 
in figure \^ then the quotient of theta functions has the following jump discontinuities on 
E 

e{uj{z) + A)\ fe{uj{z) + A) 



u{z) + B)J^ \9{u{z) + B) J 

9{uj{z) + B)J^ [9i-u{z) + B)J_' 
where A and B are arhitrary 2n — 1 vectors and Aq = Bq = 0. 
Proof. The holomorphic differentials duj are given by 

duji = ''^^} dz, 

W[Z) 

for some polynomial Pi{z) of degree less than 2n — 1 in z. This means that, under the 
action of p, dcUj becomes — dcUj. In particular, we have 



p{uj){z) = -uj{z) 



(6.10) 
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where the action of p on a; is given by (16.81) and fl6.9p . 

We first consider the jumps across the gaps Sj. Take two distinct paths from Ai to a 
point z G Sj. Assume also that both curves do not intersect S and that one extends to 
the left of S, while the other to its right. The union of these paths lifts to a loop 7 on C. 
Moreover, 7 is a linear combinations of a-cycles, i.e. 

9 

7 = ^Nitti, 

i=l 

where the A^j's are non-negative integers. 
Therefore, we have 

e{uj{z) + ^ + A) \ ^ r e{iu{z) + A) 

9{u;{z) + ^ + B)J ^ ~ \9{uj{z) + B) 

e{uj{z) + A) \ ~ 
e{uj{z) + B) ) ^"'^^ 

9 
i=l 

where Ei is the column vector with 1 in the i^^ entry and zero elsewhere. 

Now consider the jumps on the branch cuts Sj. Let ^ G Sj, then take a loop 7 on £ 
consisting of two distinct curves joining Ai to z, , both non-intersecting S; one on the left 
of the cut in Ci, the other on the right of the cut in C2. This closed loop 7 is homologic 
to the 6-cycle bj. Therefore, 

■e{u;iz) + ^ + A)\ fe{u{z) + A)\ _2.^(A, 



e{uj{z) + ^ + B) J ^ \e{uj{z) + B) 

e{-u;{z) + A) 



e{-uj{z) + B)j 

Q'fc = liki- 

This proves the lemma. □ 
We can now solve the Riemann-Hilbert problem (16. 4p . (16. 5p . Let us define 



r 

2 



2n 



Y^u{z-')-K, (6.11: 



=2 



Mz) := / dA, 
'+00 



where dA is the normalized differential of third type with simple poles at 00"*^ and residues 
±^ respectively. In addition, we write 
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Proposition 3. Let oo"^ he the points on C that projects to oo on Ci. The unique solution 
of the Riemann-Hilbert problem Ii6. 5\) is given by 



where entries of Q{z) are given by 



Quiz) = v/J^A^e-^^^) 



-10- 








z)+p{X)^- 






e [uj{z) + 1) 

(z)-/3(A)^ + 




)^(^ 


e {u;{z) - 
{z)+l3{X)-t + 


1 


e{uj{ 


e (a;(.) - I) 


1 


9 {u;{z) + f ) 



(6.12) 



621(2;) = -^z-Xie ^'—^ 7^T~7~^ 7\ ^ (6-13) 

(z) = V7^e"^(^)- 

where and ~t is a 2n — 1 dimensional vector whose last n entries are 1 and the first n — 1 
entries are 0. The branch cut of ^/ z — Xi is defined to be S/Sq. 

Proof. By using lemma [21 we see that Q{z) has the following jump discontinuities 

(©11(2;))+ = (©i2(2;))_, zeT^i, i = l,...,n 

(©12(2;))+ = (©ii(2))_, zEEi, i = l,...,n 

(©21(2;))+ = (©22(2;))_, zeT.i, i = l,...,n 

(©22(2;))+ = (©2i(2;))_, zeEi, i = l,...,n 

(©ll(^))+ - ^"^ 

(©12(^))+ - 
(©2l(^))+ - 
(©22(^))+ - 



x + 


1 


x + 


1 


X- 


1 


X- 


1 


x + 


1 


x + 


1 


X - 


1 



(©12( 


»)-, 


z 




i = n + 1, . . 


..,2n 


(©111 


»)-, 


z 




i = n + 1, . , 


..,2n 


(©22( 


»)-, 


z 




i = n + 1, . . 


..,2n 


(©2l( 




z 




i = n + 1, . . 


..,2n 



This means that Q{z) has the same jump discontinuities as in (E 



To see that Q{z) has the singularity structure given by fl6.5p . note that the function 

U+ = Q{z)e-\z), \z\ < 1 
U- = e{z)AQ-\z), \z\ > 1 

has no jump discontinuities across the branch cuts Sj. It can at only have singularities 
of order less than or equal to | at the points z^^. This means that, if it was singular 
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at zf^, then it would have jump discontinuities across Sj due to the branch point type 
singularities. Therefore it is holomorphic at the points z^^. Hence, the function Q(z) must 
have the singularity structure of the form 06.51) . 

To show that S{z) has the correct asymptotic behavior at 2; = 00, we only need to 
prove that Q{z) is invertible at z = 00. The asymptotic behavior of Q{z) is given by 



622(00) = 9(u{oo)-K-(3{X)^ + ^^e' 



Ao 
■Ao 



612(00) = 621(00) = 



where Aq = lim^^oo ^{z) - | log(2; - Ai). 

We will now show that u;(oo) = n. Let 77 be a third type differential with simple poles 
at the points Xi & C and be a holomorphic differential. Let 11* and 11* be their periods 



/ r/ = m, f 

Jai Jbi 



Now, by the Riemann bilinear relation [TB] we have 

- W+'W = 2m Res^, {t]) / f/, 

where po is an arbitrary point on C. By substituting rj = dA and fj = duj for j = 1, . . . , g, 
we see that 

Kj = - (c<jj(oo"'") — ujj{oo~)) = ujj{oo), 
where the last equality follows from (16.101) . Therefore, we obtain 

811(00) = e (/5(A)^ + e-^« 

622(00) = e(-/3(A)^ + 0e"^° (6.14) 
612(00) = 621(00) = 

Therefore Q{z) is invertible at 00 as long as 

{m^ + 1) {-m^ + D ^ 0- (6.15) 

Thus, S{z) is the unique solution of the Riemann-Hilbert problem (16. 4p . □ 
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Remark 4. In appendix F, we will show that the Wiener- Hop f factorization is solvable for 
P{X) e zM, i.e. the Riemann-Hilbert problem is solvable for these /9(A). This in turn 
implies that ^6.15\) is true for all /?(A) G iM. Define (cf. ( (5. 6j) ) 

f], = {A G M : |A| > 1 + e}. 

The function A — * /3(A) maps Vt^ onto the bounded subset A/" = {a G : < |a| < 
2^ log(2e~-'^ + 1). By continuity, the inequality l\6.15\) is valid for all a from the closure 
of M . This fact, together with the explicit formulae Ii6.12\) . Ii6. 3\) and Ii6. 1\) implies the 
uniform estimates which have been stated in ( [5. 7\ ) and used in the proof of theorem\^ 

7 The asymptotics of dlog D / and Dl(X) 

We are now ready to compute the derivative of the determinant Dl{X). First we notice 
that in virtue of fl6.1l) . equation (15.81) can be re- written as 

+^|^^ ^tr [U[iz)U-\z) {<^-\z) - as^-\z)as)]dz 
+ri(A). (7.1) 

Define 

*w-*-*w--=*-'w-3-r^(,4, "f')- 

From equations (16.31) and (I6.12p we have 

U-\z) = Aeiz)Q-\z), U'^z) = Q'iz)Q~\z)A~' + Qiz)iQ''yiz)A-\ 
where we denote A = Q{oo)A^^Q^^{oo). Furthermore, from equation (I6.2p we obtain 



1 / g \z) -i 



2 V -9-\^) 

Therefore, formula (17. ip transforms into the relation 
— log Dl{X) 



dA"" 1-A2 

dz 

+r4A). (7.2) 

We will now prove the following: 
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Theorem 6. Let s{X) be given by 

s{X) - 



a{z)dz, 



(7.3) 



z\ = l 



a{z) = tr 



dz 



where the entries of the 2x2 matrix Q{z) are given by Ii6.13\) . 
Then s(A) can be written as 



7r(l-A2)d/3 



log {e (^/3(A)^ + ^)^(/?(A)^- 



Proof. To begin with, we would like to treat a{z)dz as a 1-form on the hyperelliptic curve 
C We will show that it is, in fact, the holomorphic 1-form 



2n-l 



a{z)dz = J2 ^i^^S (O (/9(A)^ 

i=l 



0^(/3(A)^-£))da;, 



where dui are the normalized holomorphic differentials on C and di is the partial derivative 
with respect to the i*'* argument. 

Suppose this is true, then by deforming the contour of the integral (17.31) . we see that it 
can be written as 



.(A) 



TTd-A^ 



k=n -^"-k 

2n— 1 „ 2n— 1 

fc=n -^"-k j=i 
2n-l 

^ 9, log {e (/3(A)^ + e (/3(A)^ - £ 

log [e (/3(A)e^ + ^)^(/?(A)^- 



k=n 



d 
'dp 



To see that a(z)dz is given by the corresponding 1-form, let us first compute a{z)dz. We 
have 



a{z)dz = (dete(z))-M e22(^)e;i(z) -eii(z)e'22(2) 



(7.4) 



where the prime denotes the derivative with respect to z. 
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We can simplify equation fl7.4p by observing that 



eii(^ 

hi{z 
h2(z 



h,{z)9(u{z)-P{\)r-K + ^^ 
h2{z)9(u;{z)-P{\)r + ^-^) 



h2iz)e (u;{z) + + K 



-A{z) 



z — \i 



e (ujiz) + 1) 



e (uiz) - i) 



Therefore, we have 



where the ^j's are given by 

Oi 

02 

9, 



{h2{z))' {6,6', - 6,6',) , 



T 



6{uj{z)+(3{X)^ -K + - 



T 



uj{z) - (3{X)^ - K + - 



T 



uj{z) - (3{X)^ + K - - 



T 



uj{z) + (3{\)'t + ft; - - 



Now, the ^-'s are just 



2n-l 



6',dz= Y,Wi)dujk. 



k=l 



By substituting the right hand side of this equation into (17.41) we obtain 



2n-l 

a{z)dz = det Q{z)-^ ^ duk {{h,{z)fGl{z) - {h2{z)fGl{z)) 

k=l 

~- 62dk6i — 618^62 
- 6,dk6i — 6idk6z. 



Gl{z) 
Gliz) 



We would like to show that the expression 

det e{zr'{{h{z)f Gliz) -{h2{z)r Gliz)) 
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is a constant. First note that, by considering the jump and singularity structure of det Q{z), 
we have 

det 6(2) = g{z) det e {00) g {00)"^, 

where g{z) is given by (12. 3p . 

Since the 6jj(2;)'s have square root singularities at the n points z = z~^, the functions 

{h,{z)fGl{z) - {h,{z)fGl{z) 

can have at most simple poles at the points [zj)"^^, j = 1, . . . , 2n. Near each of these points, 
they behave like 

{h{z)rGl{z) - {h,{z)rGliz) = Ai + A{{z-z,y^+0{z-z,), z ^ z, 



{h,{z)fGl{z)-{h,{z)fGl{z) = Bl{z-zT^r^ + B{{z-z~'r-^+0{l), z ^ z'' 



-^3 

3 ' ' 3 

Since p(A)(2;) = — A(2), p{ijS){z) = —uj{z) and p{z — Ai) = z — \i, we have 

p{hj){z) = hliz), pie^)iz) = es{z), p{e2)iz) = e,{z) 

and 

{h{z)rGl{z) - {h,{z)fGl{z) = {h,{z)fGl{z) - p{{h,fGl){z). (7.5) 
Since the action of p on a Laurent series near a branch point \j is given by 



P 



\k=~oo / k=~oo 

by (17.51) we obtain Aq = = for all j. Hence, the function 

det e{z)-' {{h,{z)fGl{z) - {h,{z)fGl{z)) (7.6) 
does not have any pole on C To see that it does not have jumps too, let us consider 

{h,{z)fG\{z) = {h,{z)f {9M - Oid^e^) . 
The periodicity of the term inside the brackets is given by proposition O 
eM{z + aj) = 

- ^i^2(2vri(5jfc)e"2''^(2^^(^)"2''^+^^+""^ 
e2dkei{z + aj) = 6^92 

- ^2^i(27ri5j-fc)e"2''*(2^^(^)-2''^+^^+"^^), 
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where ujj{z) = J^^ duj is the j*^' component of the vector u!{z). Hence the multiphcative 
factor picked up by Gl{z) after going around a 6-cycle cancels exactly with the factor 
picked up by {hi{z))'^. It follows that the function fl7.6p does not have jumps on C too. 
Hence, they are holomorphic functions on C without any pole and must be constants. 
These constants can be computed by taking z = oo. In other words, they are given by 
(16.141) . We therefore have 

deteiz)-' {{h,{z)fGl{z) - {h,{z)fGl{z)) = a, log {o (/3(A)^ + ^) 9 (/3(A)^ - 0) 

This proves the theorem. □ 

Theorem [6l in its turn, yields our main asymptotic result. 

Theorem 7. Let be the domain of solvability ( f5. 6\} . Then the logarithmic derivative 
of Toeplitz determinant Dl{X) admits the following asymptotic representation, which is 
uniform in X & Q^. 



^logi^^A) = "Y^^+^iogK^^^^^+DK^^^^^ 



2Ji 



L^oo. (7.7) 

Here p is any real number satisfying the inequality 

1 < p < min{|Aj| : \Xj\ > 1}. 

The uniformity of the estimate (17.71) with respect to A G allows its integration over 
Q^, which yields the equation 

log (fl,(A)(l - A2)-) - Ito log (D,(s)(l - sr^) = log " '^'"^ ""I'^.P'^ " ^' 

+r(i), 

where r{L) = O (p~^) as L — oo. Taking into account (14. lip , the second term in the left 
hand side is zero. This proves Proposition [H 

8 The limiting entropy 

Observe that equation (14.101) can also be rewritten as 

SLipA)= liniT^/ e(l + e,A);^log(Z}i(A)(A2-l)-^)dA. (8.1) 
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The right hand side of this equation follows from 

lim / e(l + e,A)-5^1og(A2-l)~^dA = L lim / e(l + e,A)^^dA 
^^o+/^/,^ dA ^-*o+/r(,) 1 - A^ 



27riL lim 



'r(.) 

2X \ ( 2A 

resA=i ( e(l + e, A)- — I + resA=_i I e(l + e, A)- 



1-A2 



l-niL lim ( (2 + e) log + e log ^ 



l-A^ 



0. 



We identify the limiting entropy Si^pA) as the following double limit (cf. [I3 

lim i) 



lim 



dA. 



(8.21 



We now want to apply theorem [7] and evaluate the large L limit in the right hand side of 
this equation. To this end we need first to replace the integration along the contour r(e) 
by the integration along a subset of the set Vt^ where we can use the uniform asymptotic 
formula fl7.7p . 
Let us define 

w ■.= j-^\og{D,{\){\'-ir'^). 

The function 5(A) satisfies the following properties. 

1. (5(A) is analytic outside of the interval [—1, 1]. 

2. 5{-\) = -(5(A). 

3. 5(A) =0{X'^), A^oo. 

4. 5(A) =0(log|l-A2|), A^±l. 
Consider the identity 



/ e(l+e,A)-^log(Di(A)(A2-l)"^)dA= / e(l + e, 
Jr(e) cl^ Jr{€) 



A)5(A)dA. 



Property 1 allows us to replace the contour of integration r(e) by the large contour T' as 
depicted in figure 1, so that 



e(l + e,A)5(A)dA= 6 e(l + e, A)5(A)dA. 



Simultaneously, property 3 allows to push i? — ^> oo in the right hand side of the last formula 
and hence re-write it as the relation, 



e(l+e,A)5(A)dA 



-l-e 



5(A) 



1 + e + A A /1 + e + A 
log+ 



log_ 



1 + e + A^ 



dA 
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6{X) 



l+e 



l+e-\ / fl+e-\ 
log. 



log- 



l + e- A 



dA. 



(8.3) 



Here log^ log- ( ^^n^^ ) denote, respectively, the upper and lower boundary 

values of the functions log ( ^"^^^ ) on the real axis. We note that 



and 



log. 



log_, 



1 + e + A 



1 + e- A 



log_ ( 1 = 27rz, for all A < -1 - e, 



log_ 



1+e-A 



-27rz, for all A > 1 + e. 



Therefore, equation (18.31) becomes 
e(l + e,A)5(A)dA = 



vri 



oo 
oo 



:i + e + A)5(A)dA + TTi / (l + e- A)<5(A)dA 



2m / (l+e-A)5(A)dA, 

'l+e 



where we have also taken into account the oddness of the function 5(A), i.e. property 2. 
Recalling the definition of the function 5(A), we arrive at 

i e(l+e,A) — log(Di(A)(A2-l)-^)dA = 27ri / (1+e-A)— log (Di(A)(A2 - 1)"^) dA 

JV{e) OA dA 

(8.5) 

The estimate (17.71) can be used in the right hand side of formula (18.51) . This enables us 
to perform an explicit evaluation of the large L limit in (18. 2p so that the formula for the 
entropy S{pa) becomes 



- lim 

2 e^0+ 



^J(l + e - A)^ log(^ (/3(A)^ + e (/3(A)^ " 9)^^ 



(8.6) 



To complete the evaluation of the entropy, we need to prove the existence of this limit. 



9 Integrability at ±1. The final formula for the en- 
tropy 

We will now proof the integrability of the function 



log 
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at ±1. 

First let us denote the real and imaginary parts of the period matrix 11 by Re 11 and 
Im n. Since the Im 11 is non-singular, there exist a real vector if such that 

"e* = Imlllf 

We now can write 

i-t = (n-Ren)V. 

Let Q he a. large real number, and let nt be an integer vector such that 

Qlf = Tfi + ~^, (9.1) 

where the entries of if are between and 1. 
In particular, we have 

nt = Q (Im n)~^ "e* — if. 
Then, from the periodicity of the theta function fl6.6p . we see that 

= exp l^gV i(-e',(Imn)"^Ren(Imn)"^^) 

+ (^,(imn)-'^) 

-2nrQ [{-t , {lmU)~^ ReWf) + ("e", (ImH)"^ T'o) + {2i-t,-f)] 

x9{-{m + ^fRen + ^o + ^^U) (9.2) 

for some bounded constant "cfg- 

Note that there exists an integer vector / and real vector with entries between 
and 1 such that 

{rft + "g^)^Re 11 = T + T*. 

Therefore, we have 

e{-{rn + -ffRe n + ^0 + T^n) = e {-y + ~to+ g'^n) . 

If log^ {iQ~^ + ifo) is non-zero for all Q, then from (19. 2p we see that 



log6' {iQ~i + if 



Q'tt 



-t, (Im n)"^ Re n (Im H) ^) 



+ (Im n)-^ ^> + 2.^^^ + 0{Q-') 
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where N{Q, cq) is an integer that depends on the branch of the logarithm. It may depend 
on Q and "?o- This term arises because in the integral expression of the entropy, 

the branch of the logarithm must be chosen so that the integrand is continuous in A. We 
shall determine the asymptotic behavior of N{Q, cq) as Q — *■ 00. 

Due to theorem [HI the inequality (16.151) is true when /3(A) G zM. Therefore, we can 
apply the above result to compute the asymptotic behavior of the integrand in ()9.3p : 

" [2) 



e^,(Imn) ^ "e^) 
i ("e", (Im n)"^ Re U (Im H)"^ ^) 
2i- 



.iV(/5(A),f)+iV(/?(A),-f; 



+0(/3(A)-i ' 



(9.4) 



Since -Dl(A) in (14.111) is real and positive for A G (1, 00), and that logDi(A)(A^ — 1)"-^ has 
to be zero at A = 00 (which is needed to deform the contour to obtain (18. 6p ). we see that 
logL'i(A) has to be real for A G (1, 00). Therefore, the imaginary part of the leading order 
term in (19.41) must be zero. In particular, this means that 

/-a m-iR na m-i-\ . Ar(/3(A), j) + iV(/3(A), -j) 

(e,(Imn) Ren(Imn) e)-2 ^^^^y^ = 0(/3(A) ). 

Thus, the asymptotic behavior of the integrand in (19. 3p is 

log Hm+^) m\)-^^) ^ ^ 0{m~')) , A - l+(9.5) 

^ [2) 

The left hand side of this equation is therefore integrable at A = 1+ and we can take the 
limit e — s> in (18. 6p to obtain our final result for the entropy: 



10 Critical behavior as roots of g{z) approaches the 
unit circle 



The purpose of this section is to prove theorem [21 We shall study the critical behavior of 
the entropy of entanglement as some pairs of the roots (12.40 approach the unit circle. As 
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we discussed in section |2l in each pair one root lies inside the unit circle, while the other 
outside. In this limit the entropy becomes singular. We shall study all the possible cases 
of such degeneracy, namely the following three: 

1. the limit of two real roots approaching 1; 

2. the limit of 2r pairs of complex roots approaching the unit circle; 

3. the limit of 2r pairs of complex roots approaching the unit circle together with one 
pair of real roots approaching 1. 

When pairs of roots in fl2.4p approach the unit circle, the period matrix 11 in the 
definition of the theta function (12.101) becomes degenerate and some of its entries tend 
to zero. This will lead to a divergence in the sum (12.101) and hence a divergence in the 
entropy. It is very difficult to study such divergence directly from the sum (12.101) . In 
order to compute such limits, we need to perform modular transformations to the theta 
functions. In particular, the following theorem from [7] will be used throughout the whole 
section. 

Theorem 8. // the canonical bases of cycles {A B) and [A B) are related by 



B j \^ J V Z21 

where the matrix Z is symplectic i. e. 



hn^l 

then we have the following relations between the theta functions with different period ma- 
trices: 



Zl2^ 




Z22 ) 







— -^2ri-l 







^22 


-Zl2\ 


^21 


ZL ) 



9 

where 



(e, n) = <^ exp -meX-Zl^n + Z'Q~'Zj,^ 



9 



(e,n), (10.1) 



-Z^^fi + Z^^f)^ (10.2) 

and is a constant. The characteristics of the theta functions are related by 

e = Zj2£ + ^iV-diag(Zf2^22) 
e' = Zle + Z^,i'-diag{Zf,Z2i), 
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where (iia,g{C D^) is a column vector whose entries are the diagonal elements of CD^ . The 
new period matrix is given by 

II = {Z22II + Z21) {ZuH + Zn)'^ (10.3) 
and the normalized one forms are related by 

dn = (^{-zf^ii + zj^f^dn (10.4) 

d^F = (dcji, . . . , dcl;2n-i)^, dfi^ = {dui, . . . , du;2n-i)^, 

which is the same transformation as in U0.2\) . 

Our aim is to find a good choice of basis {A B) sucfi tliat 6{^, 11) remains finite wfiile 
some entries of 11 tend to infinity as certain pairs of roots Xj approacli tlie unit circle. Tliis 
would confine tlie divergence of tlie entropy witliin tlie exponential factor in fllO.lj) . which 
can be computed. 

10.1 The limit of two real roots approaching 1 

In this section the choice of the basis {A B) described in theorem [8] is the one shown in 
figure In the notation of theorem [H], the new basis [A B) and the old one {A B) are 




Figure 5: The choice of cycles on the hyperelliptic curve C The arrows denote the orien- 
tations of the cycles and branch cuts. 



related by 



z = 


( Z\\ Z\2 
\Z2\ Z22 ^ 




r 


-C2\ 

J 




A' = 


(ai, . . . , a2n- 


-l) 5 




=(&!,... 




A^ = 


(Oi, . . . , a2n~ 




B^ 


=(&!,... 


; ^2n-l)"^ 


(,Ci)ij = 


1, 3 > h 


(Ci; 


).. = o 


, j <i 




iC2h = 


1, , {^2)1,1- 


-1 = 


-1, 


{C2)ij = 


0, j 7^i,i 


Ci = 













;io.5) 
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The relation between the two period matrices can be found using (110.31) 

n = -CiU-^C^\ (10.6) 

To study the behavior of the entropy as the real roots \2n ^2n, we need to know the 
behavior of the period matrix 11 in this limit. Now, we have 



Wo 



A2 



4n 



lim 



A,;) 



An 



n (^-^^)- 



(10.7) 



Furthermore, as \2n ^2n ^^e integration around a„ tends the residue at ^ = 1; the 
hyperelliptic curve C becomes a singular hyperelliptic curve Cq of genus 2n — 2; the tilded 
basis of canonical cycles on this curve reduces to 



^0 



(ai, . . . , a„_l, CLn+li ■ ■ ■ J (^2n-l) ; 
(6l, . . . , bn-l, bn+1, • • • , b2n-l)'^ ■ 



(10.8) 



The holomorphic 1-forms duj tend to the following limit [T]: 



Figure 6: As \2n ^2ny integration around a„ becomes a residue integral around z = 1. 



^ Wo 

where ^Pj{X) are degree 2n — 2 polynomials determined by the normalization conditions 

dul = 6kj, j 



27riRes 



~,o 



2=i,to=-!«o(i)dci;^ 



Therefore, the 1-forms du^, k n, become the holomorphic 1-forms that are dual to the 
basis Ao on Co- Furthermore, du^ becomes a normalized meromorphic 1-form with simple 
poles at the points above ^ = 1 on £o- 

As in [T], we see that the entries of the period matrix 11 tend to the following limits: 



lim Ujk 

J- J--nr) 



n-fc, i,j^n,n 



2J2 / dcD„ 
j=i -J^ij-i 

— log \^2n - A2n| + 0(1), \2n ^ >^2n^ 
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where 11^^ is finite for z, j ^ n, n. 

Let us adopt tfie notation of theorem [8] and denote the argument of the theta function 
in the entropy fl9.6l) by that is 

e = /3(A)^±^. (10.9) 

We will now compute the behavior of the argument ^ in fllO.ip with ^ given by (110.91) . 
We have 

Lemma 3. Let ^ be given by /[lU.y\) and ^ be 

where Zij are given by ( fiO.5]) . Then as X-^n — ^ Ag^^ we have 

= /3(A)ni„±r/„ z = l,...,2n-l, (10.10) 
where r]i remains finite as X-zn Ag^^ . 

Proof. To begin with, we will need to express | in terms of the Abel map. 
Recall that the term | in (12.121) is given by 

2n 

where K is the Riemann constant. As in [6] (see also appendix D), the Riemann constant 
can be expressed as a sum of images of branch points under the Abel map. In particular, 
we have 

2n 

K = -Y,^{X2,^,). 

Therefore we have 

2n 2n 

^ = -5Z^(^7') + I]^(^2,-l) 
i=2 j=2 

Now by substituting (110. 5p into (110.21) and make use of (110. 4p and (110.60 . we see that the 
argument ^ in ^(^,11) can be expressed as follows 

(2n 2n \ 

J]cD,(^7^)-J]cDi(A2,-i) , 2 = l,...,2n-l, (10.11) 
i=2 j=i / 
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where u is the Abel map with dci; replaced by du and Ui is the z*'^ component of the map. 
We would like to show that the term 

2n 2n 

i=2 j=i 

in (110. Ill) remains finite as \2n ^2n- 

To see this, note that the set of points {zj^} must contain either one of the points 
or X2n, but not both, while {A2i-i} contains only. As — > A2„\ the terms u;„(A2n) 
and uJn{^2n) equation fllO.lip will tend to — oo. However, since they appear 

in the sum with opposite signs, these contributions cancel and the quantity 

2n 2n 

j=2 j=l 

remains finite as A2ri A2„^ 
We can therefore write ^ as 

= P{X)flin±r]i, i = l,...,2n-l 

where 77^ remains finite as A2n ~^ X2n- n 

We are now ready to apply theorem [HI to compute the theta function as A2n — ^ A^„^. 

Lemma 4. In the limit X2n Ag^ the theta function 9(^, IT) behaves like 

9{^, n) = exp (log |A2„ - A2- V'(A) + 0(1)) , (10.12) 

where ^ is given by ( liO.gj) . 

Proof. Firstly, let us use (110. ip and (110.60 to express 9{^, U) in terms of 9{^, 11), we have 

9{^,U) = ?exp vrif^n-^f ^ (^f,n) . (10.13) 

Let us now use (110. lOp to compute the asymptotic of the exponential term in (110.130 . We 
obtain 



The behavior of the entries in 11 ^ can be calculated by computing the determinant and 
the minors. We have 

(n^')^. = 0(1), A2n^A2„\ i,j^n 

(n-i) = 0(log-i|A2„-A2-^|), A2„^A2-„\ j^n 

\ / nj 

(n~') = 7rzl0g-MA2n-A2„l|+0(l0g-'|A2n-A2„l|), X2n ^ X^^ . 

\ / nri 
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Therefore, equation (110.141) becomes 

Y.i^~'),M3 = l0g|A2n - \2n\P'W + 0(1), \2n ^ A^'^^ (10.14) 



Next, we will use the definition (12.101) of the theta function to compute its limit as A 



2n 



^2n- We have, 



9{^,Il) = ^ exp vrz ^ n^fcrn^mfc + 27rz^ (^/5(A)nj>, ± r/j 

+27rin„„ {ml + 2/3(A)m„) ± 2r/„m„ . (10.15) 

Since 

lim Re(27rzn„„) = — oo 

and /3(A) is purely imaginary, we see that in the limit only the terms with m„ = contribute 
to the sum. Therefore, equation (110. 15p reduces to 

lim e{iu) = eU',u') (10.16) 

= (^1, • • • 5 ^ni • • • 5 ^2n-l)'^, 

where the ^„ in the above equation means that the n*'* entry of the vector is removed. The 
period matrix 11° is an (2n — 2) x (2n — 2) matrix obtained by removing the n^^ row and 

n*'^ column of the period matrix fl. Thus, the theta function 9 H'' j remains finite as 
A2n ^2n- This fact, together with (110.141) . shows that 6'(^,n) behaves like 

e{^, n) = ^exp (log \X2n - X2n\P'W + 0(1)) 6 (l", 6°) , X2n ^ A^"^^ 

Since 9 ^,^",11''^ and ^ remain finite as X2n A^„^, the above equation becomes (110.121) . 
This proves the lemma. □ 

Finally, by substituting (110.121) into (19. 6p . we have 

, 1 Ti ^(/^(A)^ + i)g(/3(A)^-i) ,, 
= 2/ log dA 



(log|A2„-A2-„V'(A) + 0(l))dA 



Since 



/3^(A)dA = --, 
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A 



2(j+l) 



A 



2n 



A 



2n-l 



A 



2n+l 



A 



2n+2 



A 



2(2n-j)-l 



A 



2(2n-j) 



Figure 7: Two pairs of roots, labelled according to the ordering (12.51) . approaching the unit 



circle in the critical limit. We have A2(j+i) 
respectively. 



A 



2(2n- 



A 



2n 



A 



2n+l 



and X2j+i A2(2n-j) 



we arrive at the following expression for the entropy of entanglement 



S{Pa) = -- log |A2n - A2"„'| + 0(1), A2n ^ ■ 



10.2 The limit of complex roots approaching the unit circle 

We will now study the case when 2r pairs of complex roots approach each other towards 
the unit circle. Let A2j+i be a complex root with n — r<j<n — 1. As we discussed in 
section[21 A2j+i, l/A2j+i and l/A2j+i are roots too. The ordering (12. 5p implies (see figure[7]) 



A 



2(i+i) 



A 



2j+l 



A 



2(2n-j) 



1/A 



2(2n-i)-l 



A 



2(2n-i) 



2(j+l) ^2(2n-i)-l 



1/A 



2j+l- 



(10.17) 



The critical limit occurs as A2(j+i) 



A 



2(2n-j)-l- 



From the relations (110.171) this implies 



2i+l 
2(2n-i)-l- 



A 



2{2n-i)- 



Thus, in what follows we shall mainly discuss the limit A2(j+i) 



10.2.1 Case 1: r < n 

We now choose the tilded canonical basis of the cycles {A B) as in figure [HI Namely, we 
have 



= ttj, j < n — r, j > n + r — 1 
= bj, j < n — r, j > n + r — 1 
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^n—k 



-'n—k 



M+k 



n+k-2 

bfi-k — bn+k-1 + , k = 1, . . . ,r 

j=n—k+l 
n+k 

K+k - K-k-i + ^ Oj, , A; = 0, . . . , r - 1 

j=n—k—l 
n+k— 2 n—k—1 

bn-k - 5^ a, - ^ (a^. _ 26^.) , ,k = l, 

j=n~k j=n—r 

bn+k + Yl (-l)""'^'' -'2b,), k = 0,...,r-l. 

j=n—r 



(10.18) 




Figure 8: The choice of cycles on the hyperelhptic curve C The arrows denote the orien- 
tations of the cycles and branch cuts. 

We will show in appendix E that this is indeed a canonical basis of cycles. We can 
partition this basis as follows: 



A' 



af 
a- 



ttj, l<k<n — r — 1 



"n—r+j — li 



1, 1 <k<2r 



ttj, n + r<k<2n — 1. 



The relations among the 6-cycles and the untilded basis are analogous. 
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If we write this relation in matrix form as in theorem [HI then the corresponding trans- 
formation matrix is given by 



Zn Zi2 
Z21 Z22 



where the blocks Zij can be written as 



Zii 








Cij 








^ij-^n—r 



where In-r-i is the identity matrix of dimension n — r — 1 and the Cj/s are the following 
2r X 2r matrices: 



(Cii),, 
(^22)^ 



1 k + 1 <l <2r -k, 

otherwise, 

1 k<l<2r-k + l, 
otherwise, 

hi - Si,2r-k+i 1 <k,l <2r 

1</<A;-1; 
-1 k <l <2r -k, 

2r-k + l<l 

l<l<2r-k, 
otherwise, 

2(_l)fc-« 1</<A;-1 
otherwise, 
hi - 2 iC2i)ki , r + l<k<2r 



6,. 



kl 



l<k<r 
r + 1 <k <2r 



, l<k<r 

r + l<k <2r 
1< k<r 



These are matrices of the form 

/O 1 1 
1 



C 



12 





... 



Vi 1 1 
/o 

... 1 



1 



1 



1 



1 



1 0\ 






... 

1 1/ 



1\ 





(10.19) 



Vl ... 0/ 



48 



A. R. Its, F. Mezzadh and M. Y. Mo 



/-I 
1 

-1 



a 



21 



1 

\ 

/ 1 



a 



22 



V 



1 





1 

2 
2 






-1 
-1 
-1 

-1 
-1 







-2 
-2 






-1 -1 Q\ 
-10 























0/ 








1 





Q\ 





... 








1 







1/ 



As in section [TOl some holomorphic 1-forms dujj will become meromorphic as the roots 
approach the unit circle. 

In this case, the holomorphic 1-form dcDj, n — r < k < n + r — 1 becomes a meromorphic 
1-form with a simple pole at A2(j+i). AH the other holomorphic 1-forms become normalized 
holomorphic 1-forms in the resulting surface. 

In particular, we have the following: 

Lemma 5. The entries of the period matrix U behave like 



^20 + 1) 


lim 


-j)- 


-1 


^2(j+l) 


lim 


-j)- 


1 


^2(j + l) 


lim 

-^-^2(2n- 


-j)- 


1 



4, ^y^J 
n°^-, j > n + r 



1, j < n — r 
n — r<j<n + r — 1 



■ log A2 



(i+i) 



A 



2{2n-i)-l| , 



(10.20) 



where 11^^- are finite. 

Let us now consider the behavior of the terms in (110.11) . 
Lemma 6. Let C, be given by 1^10. 9\} and ^ be 



T \T 
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where Zij are given by 1^10. 19\) . Then in the limit X2{j+i) ^2{2n~j)~i we have 

= r]f, i > n + r — 1, i < n — r 

ii = eip{\)-fi + rit, n-r <i<n + r-l, (10.21) 

ej = 1, i < n 

€{ = —1. i > n 

where rjf remains finite as A2(j+i) ^2{2n-j)-i- 
Proof. Let 

/O 0\ 

Zf^U - Zl^ = {hr - J2r)Dr o] +W 

\0 0/ 

Dr = diag(7n-r, ln-r+1, Jn-r+l, Jn-r) , (10.22) 

where is a matrix that remains finite as A2(j+i) \2(2n-j)-i ■ Then from (110.21) and 
(110.41) . we see that ^ is given by 

n-l 

ii = /3(A) Wn+j,i ± ^, i>n + r-l, i<n-r 

n-l 

ii = ei(3iXhi + f3iX)J2Wn+j,i±^, n-r<t<n + r-l, (10.23) 



where 



2 



€i = 1, i < n 
Ci = —1, i > n 

2n 2n 



Let A2{j+i), A2{2n-j)-i and A2j+i, A2{2n-j)5 n — r<j<n — Ihe the pairs of points that 
approach each other. From their ordering we have A2{j+i) = X2{2n-j) -^sj+i = A2(2n_j)_i- 
For each fixed j, the point A2j+i is a pole of doj2n-j-i, while A2{2n-j)-i is a pole of dojj. 
Therefore, the Riemann constant behaves like 

2n ^ 

^u;i(A2j-i) = -7i + 0(1), A2(j+i) ^ A2(2n-j)-i n-r<i<n + r-l. 
i=i 

Moreover, among these 4 points there are exactly two points of the form z^^ for some k. 
However, since Zk are the roots of a polynomial with real coefficients, if \j = z^^ for some 
fc, then its complex conjugate \j is also of the form z^,^ for some k' . This means that either 
of the following is true: 
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1. Both A2(j+i) and \2j+i are of the form ^, 

2. Both \-2{2n-j) and A2(2n-j)-i are of the form z'^^, 
Either way, we have 

2n 



E 



-7i + 0(1), n — r<i<n-\-r — 1. 



Therefore, we can rewrite (110.230 as 



Vi 



t > n + r 



I < n — r 



eiP{X)'ji + r]f, n — r<i<n + r 



where rjf" remains finite as A2 



(i+i) 



A 



2(2n-j)-l- 



□ 



We now compute the behavior of the theta function 6'(^, 11) in this hmit. 

A2(2n-j)-i; u — T < j < 71 + T — 1, the theta function 



Lemma 7. In the limit A2(j+i) 
9{^,I\) behaves like 



^(e,n)=exp 27r^/?^(A) ^ 7i + 0(1) 



(10.24) 



j=n~r 



where ^ is given by ( liO.gj) and 7j by filO. 
Proof. From (110. ip we see that 



e{^,U) = ^exp[me {Zl,U-Z^,) ] \ (C, H) 



7T , 



e 
Ve'l 



(10.25) 



where the characteristics on the right hand side are obtained by solving the hnear equations 

diag [ZI2Z22) 
diag (^n^2i) 

The solution of this system is 



^21^ ~^ Zj^E . 



mod 2, j = l,...,2n- 1 

1 mod 2, n — r<j<n — 1] 
mod 2, otherwise. 



(10.26) 



Note that, from (12.111) and the periodicity properties of the theta function proposition [21 
characteristics that differ by an even integer vector give the same theta function. That is 



e 



(e, n) = 9 



e' + 2li 



2n-l 
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We will now compute the exponential term of (110.251) . By performing rows and columns 
operations on — Z22, we can transform its determinant into the form 

I l^n-r-l 

det f Zfafi - Z^A = det SD^ \ +W' 

VV 0„_ 

0, 1 < ? < r, 

r + 1 < i < 2r, 



for some matrix W that remains finite as \2{j+i) — ^ A2(2n-i)-i- 

This means that the leading order term of the determinant is of the order of 7^. 
That is 



(i+l) ~^ ^2{2n~j)-l, 



n-1 

Vr = W H 7fc, 
k=n—r 



where the notation 0(7[ ) means 

Oiir') = O (U^r-^ ' ^a. < r - 1, (10.27) 

Furthermore, W" is the determinant of the (2n — r — 1) x (2n — r — 1) matrix formed by 
removing the (n — r)*'^ up to the (n — l)*'^ rows and columns in W. 

Similarly, we see that the minors of — Z22 cannot contain more than r factors of 

7. In particular, this means that the inverse matrix — ^'^g j is finite as A2{j+i) — > 

A2{2n-i)-l- 

Therefore the inverse matrix (^Z]^^ — ^22 j behaves like 

(^Z^2^-Zj2y' = X' + X~' + Oi^-'), A20+I) ^ A2(2n-,)-l, (10.28) 

where is a term of order —1 in 7j and X° is a finite matrix. 
From (110. 28p and (110. 22p . we see that the leading order term of 

(zf^n - (zf^n - = i2n-i (10.29) 



gives the following 





X° I ihr - J2r)Dr | = 
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while the leading order term of 

Z^n — ) ( -Z^^n — ) = hn-i 



gives 





{hr - J2r)Dr | X° = 0. 





This implies that 



K = Xt2n-j-i. l<^<2n-l, n-r<j<n + r-l 
= ^2n-i-ij. n-r<i<n + r-l, l<j<2n-l. 

The leading order term of the bilinear product in (110.251) then becomes 

-1 /O 

\0 ' 

+0(1), ^2{j+l) A2(2n-i)-l, 

ej = 0, i < n — r, i>n + r — l, 

= 1, n — r < i < n, 
Si = —1, n<i<n + r — l 

(K-r^l 

= Dr 

V o„„ 

Let us denote V by 

/O 

p = X'M (/2. - J2.) 

\0 

Then constant term of (110.29^ gives the following 

/O 0\ 

(/2. - J2r) D„ + X°W^ = hn-l. 

\0 0/ 

By applying (110.301) to the above, we see that the entries of V are related by 
'Pi,] ='P2n-~i~i,j + + ^2n~i-i,j, u-T <l <n-l, u-r < j <n + r 
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By substituting this back into fll0.3ip . we see that the the exponential factor in (110. 25p 
behaves hke 



n-l 



(Zf,ll - Zl) ZU = 2/3'(A) 7, + 0(1). 



^10.32) 



j=n—r 



We will now show that the limit of the theta function with characteristics remains finite. 
By using the definition (12.101) . we have 



9 



(e,n) = Yl 



exp 



e 
~2 



n-l 



+4 



j=n—r 
m2n-j-l + 



e 



' ■ 2 

^2n-j-l 



X - 2/5(A) + m2„_j_i + 



0(1) 



As before, since (3{\) is purely imaginary, only terms such that 



m» + 



m2n-j-l + 



0, n — r<j<n — 1, 



contribute. Recall that from (110. 26p we have Sj = e2n-j-\ = 0, therefore 

mj = m2n-j-i =0, n — r<j<n — 1. 

Thus, as before, the theta function with characteristics reduces to a 2n — 2r — 1 dimensional 
theta function 



lim 



e 



-^2U + l)^-^2{2n-j)-l 



(e,n) = e(eo,n°^ 



[10.33) 



where the arguments on the right hand side are obtained from removing the [n — r)*^ up 
to the (n + r — 1)*'^ entries and that 9{^^, IT") is finite in the limit. 

By combining (110. 32p and (110. 33p . we see that the theta function Oi^^Tl) behaves like 



n-l 



^(e,n) = ^exp 27rz/52(A) ^ + 0{l)\ e{i\W) 



j=n-r 

This concludes the proof of the lemma. 

Finally, from lemma [7] we see that the entropy (19.60 behaves like 

n-l 

3 



□ 



n— 1 

S{Pa) = ~q log l-^ZO'+l) - A2(2n-j)-l| + 0(1), A2(j+l) A2(2n-j)-l- 



j=n—r 
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10.2.2 Case 2: r=n 

We will now consider the case when r = n. That is, all roots are complex and they all 
approach each other pairwise. The canonical basis will be chosen as in (110.181) but with 
r = n — 1, (not n) while the last elements in the basis are given by 

In other words, we have 

n+k~2 

(in-k = bn-k — bn+k-1 + , k = 1, . . . ,12 — 1 

j=n—k+l 
n+k 

a-n+k = hn+k-K-k-i+ ^ Oj, ,/c = 0,...,n-2 (10.34) 

j=n—k—l 
n+k— 2 n—k—l 

bn-k = bn-k - ^ (^j- (-1)"^^"^ («i - 2&j) , , = 1, . . . , n - 1 

j=n-k j=l 
n-k-2 

bn+k = bn+k + {a, -2b,), k = 0,...,n-2 

0'2n-l = b2n-l, &2n-l = ~0,2n-l- (10.35) 

As before, we can partition the basis as follows: 



A 



7.1 1 



a] = dj, l<k<n-r-l (10.36) 

0'2n-l- 



Furthermore, the 6-cycles and the untilded basis are connected by analogous relations. 
In the notation of theorem [8] we have 

i) - 4t)-iz" z') it 

ij J V-°/ \^21 ^22 J \D 

where the transformation matrix Z can be written in block form according to the partition 
( 110361) : 

where Cij are 2(n — 2) x 2(n — 2) matrices defined as in (110. 19p . and E is given by 

Ei^ = 0, i = 3 
E\2 = 1, E21 = —1. 
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By deformation of the contours, we see that the cycles dj become close loops around A2(j+i) 
in the critical limit. 

Let do be the closed curve that becomes a loop around A2 as A2 \4n~1 (see figure [9]). 
We have 

2n-2 

do = —b2n~l + Clj 
J=l 

n-1 2n-2 

In particular, this means that in the limit, the 1-form ibj will have a simple pole at A2(j+i) 




Figure 9: The curve going around A2. 

with residue and a simple pole at A2 with residue (— 1)'''"^^2^ for 1 < j < n — 1, {—^y^^ 
for n < j <2n — 2 and — for j = 2n — \. Thus, we arrive at the following 

Lemma 8. The entries of the period matrix behave like 

lim liij = n°., i^j, i,i^2n-l 
^ lim % = 7,+n;,, l<3<2n-2 

lim n2n-l,2n-l = '^l2n-l + '^ln-l,2n-l 

-^2(j + l)^-^2(2„_j)-l 
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lim n,-2„,„i = (-iy72n-i + n;2„_i, l<J<n-l 

lim n,-2„,_i = (-iy+S2n-i + n°2„_i, n<j<2n-2 

-'^2{j + l)-»-'»2{2n-j)-l 

n2n-lj = nj,2n-l 

Ij = — log I -^SO'+l) - A2(2n-i)-l I , 

where tl^j are finite in the limit A2(j+i) — > A2(2r!.-j)-i- 

In this case, the argument in fllO.ll) behaves as follows. 
Lemma 9. Let ^ be given by UU.9\) and ^ be 

where Zij are given by 1^10.31^ . Then in the limit A2(j+i) \2{2n-j)-i we have 

li = aiPiXhi + r/f , l<t<2n- 1, (10.38) 
ai = (1 + (-1)^+1), l<i<n-l 
ai = -(1 + (-1)^+^). n<i<2n-l 

where rjf remains finite as A2(j+i) A2(2n-j)-i- 

Proof. In this case the matrix 2^^n — takes the form 

Dn^i = diag(7i,72, . . . ,72,7i) (10.39) 

Dn-i = (-71,72, ... ,72, -7i), 

where W is a finite matrix as A2(j+i) — > A2(2n-j)-i- 
Therefore, ^ behaves like 

n-l 

i=i 

(Ti = (l + (-l)'+^), l<i<n-l 
a, = -(1 + (-1)^+^), n<i<2n-l 

2n 2n 
j=l j=l 

As in section 110.2. H the leading order terms of ^ are zero. We can therefore rewrite ^ as 
where r]f are finite in the limit. □ 
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The behavior of the theta function for this case is given by 

Lemma 10. In the limit )^2{j+i) ^2{2n-j)-i, 1 < J < 2n — 1, the theta function 6{^,Il) 
behaves like 

9{^, n) = exp (^2mP\X)J2^, + 0(1) j , (10.40) 

where ^ is given by UU.9\) and by lemma [3 
Proof. As in section 110.2.11 from (110.11) we have, 

^(e, n) = ^exp (^t~e {zlii - z^Y^ e [^^,] ii, n), (io.4i) 

where the characteristics on the right hand side are given by the same formula as before, 
with r replaced by n — 1: 

Ej = mod 2, j = 1, . . . , 2n — 1 



-J 

, _ ( I mod 2, 1 < j < n - 1; 
^ I mod 2, otherwise. 



Since there is no non-zero matrix Xq that is independent of 7^ such that the leading order 
term of 

is zero, we can write the inverse matrix ^Z^H — Z22^ as 

^Z^n - Zjgj = X~^ + 0(7,"^), A2(j+1) A2(2n-i)-l- 

where X~^ is a term that is of order —1 in the 7^. 

Then, the leading order term of the bilinear product in (110. 4ip is 

(Zf^n - Zj,) ZU = m^^D^X-' ( ^^'"-^ ^^"-^^ \ ) Dna 

+0(1), A2(j+l) A2{2n-j)-l, 

a, = (1 + (-1)'+^), l<i<n-l 

ai = -(1 + (-1)'+!), n<i<2n-l (10.42) 
Dn = diag(7i,72, ... ,72,71, 272n-i)- 
Let be the leading order term of 11: 

/ 



27 



2n-l 
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Equation (110.421) can now be rewritten as 

+0(1), A2(j+1) A2{2n-j)-l 

Ci = 1, 1 < i < n - 1 (10.43) 

ej = —1, n < i < 2n — 1 

The constant term of 

^Z^n — j (^^^n — Z22^ = hn-i 

now gives 

^ [ I J 

By substituting this back into (110.431) . we obtain 

_1 2n-l 

TT^f^ (^1^2^ - ^^2) ^12^ = X] log I -^aO+l) - A2(2n-i)-l | + 0(1). 



To complete the proof, note that in this case, the theta function in the right hand side 
of fllO.411) becomes 1: 



-^2(i+l)^^2(2n-j)-l 



hm 6 

^■2(3 + 1) 

Therefore, we have 



(e,n) = i. 



(2n-l \ 
X] + ) ' -^20+1) A2{2„-i)-l. 

This completes the proof of the lemma. □ 
Finally, by substituting (110.401) into (19. 6p . we find that the entropy behaves like 

^ 2n-l 

S{Pa) = -nYl log l^2{i + l) - A2(2n-i)-l| + 0(1), A2(j + 1) 



1) A2(2n-jVl- 
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10.3 Pairs of complex roots approaching the unit circle together 
with one pair of real roots approaching 1 

The canonical basis used in this section is shown in figure [10 

Sfc = —bk + bk-i, k < n — r, k > n + r bo = 

2n~l n+r~l 

bk = ttj — ttj, k < n — r 

j=k j=n—r 

n+k~2 

O-n-k = bn-k — bn+k-l + 0,j , k = 1, . . . ,r 

j=n—k+l 
n+k 

dn+k = bn+k - bn-k-1 + ^ k = 0, . . . ,r - 1 

j=n—k—l 

2n-l n+k-2 n-k-1 

bn-k = bn-k + (-1)'^' 5Z % - 5Z % - H {-ly-^-' (a, - 2b,) , k = l,...,r 

j=n+r j=n—k j=n—r 

2n-l n-k-2 

K+k = bn+k + (-1)'^-' XI + {-ly-^-^ (a, - 26,) , = 0, . . . , r - 1 

j=n+r j=n—r 
r-1 

O'n+r = bn-r-1 ~ bn+r + ^ ^ {'^bn+j + Ctji+j ~ '^bn-j-l + 

i=o 

2n-l 

bk = X^ ctj , k > n + r. 

j=k 

In the notation of theorem [8], the two bases are related by 

-D / V-°/ V^21 ^22/ V-D 

0' 

11 = I Cn 



-CT'"^ 

I (10.44) 

y31 y32 _(-7n.-r-l 



^12 



'2 

CT'"' Wi3 

^21 = I C21 W23 

c^'^-'-i 

0~ 

J22 = 1 C22 
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Figure 10: The choice of cycles on the hypereUiptic curve C The arrows denote the 
orientations of the cycles and branch cuts. 



where Cij are defined in fll0.19p and are k x k matrix with entries defined as in fllO.51) . 
All the entries of the matrices U^^ are 1, while the entries of V^^, V^^ and lA^^ are defined 



m 



- A.,|'-1V+^ -7-32 
■'■ij ~ ^) ' ^i,2r-j+l ~ ■'■ij ' 



■ «J 

;32 



^il^j,n— r— 1 

2(-iy5a 



Performing the same analysis as in section 110.2.11 we arrive at 
Lemma 11. The entries of the period matrix IT behave like 



lim 

-^2(j + l)^-^2(2n-j)-l 

lim 

^2{j + l)^^2{2n-i)~l 



lim 



-^2{j+l)^>'2{2n-j)-l 



7i 



4, 



J > n + r, J < n ~ r 
7j + n°^-, n-r<j<n + r 



log 1^2(^+1) A2(2n-i)- 



where n^^- are finite. 

In this case, the argument ^ is given by the following 



^10.45) 
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Lemma 12. Let ^ be given by I^10.9\) and C, be 



^12^ + ^22)"^ 



where are given by Then in the limit A2(j+i) — > A2(2n-j)-i we have 



i > n + r, i < n — r 



6 = Vi 

= ei/3(A)7i + r/f , n-r <i <n + r -1 



(10.46) 



1, i <n, ei 



2 > n — 1, 



where r]^ remains finite as A2(j+i) — > A2{2n-j)-i) n — r<j<n + r. 

The proof of this lemma follows from exactly the same type of argument as in section 
ll().2.1[ 

We will now compute the limit of the theta function. 

Lemma 13. In the limit X2{j+i) X2{2n-j)-i, n — r < j <n + r, the theta function 6{^,Il) 
behaves like 



n-l 



e{^, n) = exp I 2m(3\X) ^ 7, + f3\X)ln+r + 0(1^ 



^10.47) 



where ^ is given by U0.9\) and 7^- by ^lO.^^ - 



Proof. The characteristics in the theta function in (110. ip are once more given by (110.261) . 



The matrix — can now be written as 



/0„_._i \ 

{l2r-J2r)Dr 

7n+. 

V On-r~lJ 

d.iag(7n— r; Tn— r+l; • • • ; Tn— r+1; Ti— J")' 



+ w 



where W is finite in the limit and On-r-i is the zero matrix of dimension n — r — 1. 

As in section 110.2.11 by performing rows and columns operations on the matrix ^^j^II 
Z22, we see that the determinants has the following asymptotic behavior: 



det I Zfafi - Z^2 



In+rT^r + 0{^l), A2(j+1) ^ A2(2r!,-j)- 
n-1 

w n 

k=n—r 



where the notation 0(7[) was defined in equation (110.271) and W" is some constant. 
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The inverse matrix ( — ) can now be written as in fll0.28p : 



2(2n-j)-l, 



where the entries of the 2r dimensional matrix satisfy (110.301) with {X^)n+r,n+r = 0, 
and is a matrix of order —1 in the 7j with {X~^)n+r,n+r = In+r- 

Following exactly the same analysis in section 110.21 we see that the leading order term 
in the exponential factor in (110.251) is 

(Z?2n - Zj,^ Zli = P\X) I 2 J2 7. + 7n+r J + 0(1). 

\ j=n-r / 

We now look at the term 9 {i^^ (110.11) . As in section [10.2.11 we see that the theta 
function becomes 2n — 2r — 2 dimensional: 



lim 6 

A2{j+l)^-^2{2n-j)-l 



(e,n) = e(e°,n°^ 



where the arguments on the right hand side are obtained from removing the (n — r)*'* up 
to the (n + r — 1)*'* entries. 

Therefore the theta function 6 (^, 11) behaves like 

^(e,n) =^exp (^2mP\\) (^2 7, + 7n+. j + 0(1) j ^(e°,n°). 

This completes the proof of the lemma. □ 
By substituting (I10.47P into (19.61) . we see that the entropy is asymptotic to 

1 ""^ 1 

S{Pa) = -3 X] log 1^20+1) - A2(2n-j)-l| - g log |A2 
j=n—r 

+0(1), ^2(j+l) ^2{2n-j)-l- 

This concludes the proof of theorem [21 



Appendix A. The density matrix of a subchain 

Let be a basis of the Hilbert space of a system composed of two parts, A and B, 

so that Ti. = TYa ® '^b- The density matrix of a statistical ensemble expressed in the basis 
is a positive Hermitian matrix given by 

PAB = ^Cjk \iJj){^k\ , 
jk 
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with the condition tiABPAB = 1- Let us introduce the operators S{j, k) and S{j, k) defined 
by the relations 



SU, k) 



S{j,k)S{kJ) = 6ji \tlJj){tlJj\ and S{j,k)S{kJ) = 6ji 
(In this formula repeated indices are not summed over.) Clearly, we have 

Cjk = trAB [pab S{kJ)] . 



Let us now suppose that the Hamiltonian of our physical system is (13.101) and that 
the subsystem P is composed of the first L oscillators. Then a set of operators S{j, k) 
for the subchain P can be generated by products of the type 11^=1 where Gj can be 
any of the operators {cj, c], cjc^, c^c]} and the CjS are Fermi operators that span Ha', it is 



straightforward to check that S{k,j) = ( 11^=1 ^i) • then have 



Pa = trp 

All the S{l,k) 

= E 

All the Sil,k) 
All the S{l,k) 



L \ f 



^ t 

tiB (pab) I n^'i 
L \ t 

PAB ( Y[Gj 



J = l 



Since Pab = |*g)(*g|5 this expression simply reduces to 



PA= E (*slfn^^) \^.)tlGr 

Allthc5(«,fc) \i=l / j=l 



The correlation functions in the above sum can be computed using Wick's theorem (13.91) . 
Finally, if the correlations of the cjs are given by (14. 5 p and (14. 6p . we immediately obtain 
formula (14. 7p . 



Appendix B. The correlation matrix Cm 

The purpose of this appendix is to provide an explicit derivation of the expectation values 

(*g|mjmfc|*g) (B.l) 
when the dynamics is determined by the Hamiltonian (13.101) . 
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First, we need to diagonalize Ha, which is achieved by finding a hnear transformation 
of the operators hj of the form 



M-l 



i=o 

such that the Hamiltonian (13.101) becomes 

A/-1 



Ha=Y, lAfcl vW + C, (B.3) 

k=0 

where the coefficients gkj and hkj are real, the r/^s are Fermi operators and C is a constant. 
The quadratic form (13.101) can be transformed into ( IB. 31) by (IB. 21) if the system of equations 

[%,i7j-|Afc|r/fc = 0, A; = 0,...,M-1 (B.4) 

has a solution. Substituting (13.101) and (IB.2p into (IB.4p we obtain the eigenvalue equations 



\^k\gkj = {gkiAij — hkiBij) , 

|Afc|/ifcj = ^ {gkiBij - hkiAij) , (B.5) 
where A = aA — 21 and -B = 0:7-8. These equations can be simplified by setting 

= 9kj + hkj 

i^kj = gkj - hkj, (B.6) 
in terms of which the equations (IB.SD become 

(A + 5)0, = |A,|V, (B.7) 
(A-i?)V, = |A,|(/>,. (B.8) 

Combining these two expressions, we obtain 



(A-5)(A + B)0,, =|Afc|'0, (B.9) 



(A + 5)(A-i?)V>, =|Afc|'V.. (B.IO) 

When Afc 7^ 0, (pk and |Afc| can be determined by solving the eigenvalue equation ( IB. 91) . 
then 1/?^ can be computed using ( IB. 71) . Alternatively, one can solve equation ( IB. 101) and 
then obtain 0^ from ( IB. 81) . When A^ = 0, (pk and t/>^ differ at most by a sign and can be 
deduced directly either from (1B.7P and (IB.Sp or from (IB.QD and (IB.lOp . 
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Since A and B are real, the matrices {A — B){A + B) and {A + B){A — B) are symmetric 
and positive, which guarantees that all of their eigenvalues are positive. Furthermore, the 
0^s and V'fcS can be chosen to be real and orthonormal. As a consequence the coefficients 
Qkj and hkj obey the constraints 



A/-1 



^ {gkjgki + hkjhki) = Sji, (B.ii) 

fc=0 
M-l 

{gkjhki + hkjgki) = 0, (B.12) 



k=0 

which are necessary and sufficient conditions for the r/^s to be Fermi operators. 

The constant in equation (IB. 31) can be computed by taking the trace of Ha using the 
two expressions (13.101) and flB.3p : 

M-l M-l 

fc=0 fc=0 

Therefore, we have 



C = -5^Mfcfc-2-|A,|) 



k=0 

We are now in a position to compute the contraction pair fIB.ip . Substituting flB.6p 
into flB.2p we have 

^fc = - E {(f)kjm2j+i - i^pkjm2j) . (B.13) 

j=0 

Since the 0^'s and i/^kS are two sets of real and orthogonal vectors, (1B.13P can be inverted 
to give 

M-l 

m2j = i 



k=0 
M-l 

:j+i = Yl '^^^ ip'' + ^fe) • (^-^^^ 



m2_ 

fc=0 



Since the vacuum state of the operators rjk coincides with |^'g), the expectation values (IB.ip 
are easily computed from the expressions (lB.14p and (lB.15p . We have 



M-l 



(*g| maj-msfc |*g) = V'ZjV'fc = "^ifc, (B.16) 

1=0 
M-l 

(*g| m2j+im2k+i |*g) = Y ^^^^^^ ^ 



1=0 
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and 



M-l 

(*g| m2jm2k+i |*g) = i ^ ^ij(l)ik, 

1=0 

(*g| m2j+im2k |*g) = i^ik(pij- 

1=0 



(B.18) 
(B.19) 



Finally, by introducing the real M x M matrix 



M-l 



lj(Plk-, 



j,k = 0,...,M-l 



1=0 



and combining the expressions ( 1B.16|) . (jBTTZj), ( IB. 181) and ( IB. 191) we obtain 

(*g| mjiTik |*g) = Sjk + i{CM)jk, 

where the matrix Cm has the block structure 

^ Cu Ci2 ■ ■ ■ CiM ^ 
ri _ ^21 C22 ■ ■ ■ C2M 

— 

\ Cmi Cm2 ■ ■ ■ Cmm ) 



(B.20) 



:B.2r 



(B.22) 



with 



(B.23) 



iTM)jk 

-{TM)kj 0. 

We call Cm the correlation matrix. It is worth noting that because of the definition (IB.20p . 
the matrix Tm contains all of the physical information relating to the ground state of Ha- 



Appendix C. Thermodynamic limit of the correlation 
matrix Cm 

In this appendix we prove the following 

Lemma 14. Let he the Hamiltonian l{3.10\) and consider the correlation matrix ^B. 22\) 

(C.l) 

where Too[$] is the semi-infinite block- Toeplitz matrix with symbol 

$ = 

where the function g{z) is defined in Ii3.15\) . 



associated to Ha- We have 



lim Cm = Too[$], 

M— »oo 



g (e*^) 

-9-' {e'') 
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Proof. From the definitions f l3.15p and fl3.16p we have that 



9 {e-n = 9 (e^^) = 9'' {e^ ■ 
Thus, from equation ( IB. 231) it suffices to show that 



lim {TM)jk = 1^ r 9 (e*') e-'^^-'^'dO, (C.2) 



where g{z) is defined in (13.151) . 

The first step consists in determining the vectors 0^ and i/'fc; and the numbers via 
the eigenvalue equations (IB. 71) . ( IB. 81) . ( IB. 91) and (IB. 101) . If we use the definitions ( 13.121) . 
we can write 

(A + B)jk = a{j -k)+-fb{j -k) and (A -B)jk = a{j - k) - -fb{j - k). 

Two arbitrary circulant matrices commute and a common set of normalised eigenvectors 
is given by 

z^fc,- = -^i^=|^, j,k = 0,...,M-l, (C.3) 
VM 



where the index j labels the component of the k-th eigenvector. As a consequence, the ipf^ 
are a set of common eigenvectors of both {A + B){A — B) and {A — B). Now, combining 
equations ( IB.SP and ( IB.lOp we can write 



M-l 



- - 7Ki - 0] i^ki = ^k^kj = \K\ <Pkj, (C.4) 

1=0 

with 0^ = i/j^A'^/ |Afc|. Because both 4)j^ and i/j^ are normalized, A'^/|Afc| must be a 
complex number with modulo one and we can set A^ = A^. The eigenvalues A^ can be 
computed by directly substituting the eigenvectors ( 1C.3P into the left-hand side of ( 1C.4P 
and using the parity properties of the functions a{j) and b{j). We obtain 



Ef=:llLy2 - 7Kj)) e''^' if M is odd 

E -I'm/s-i - iKj)) e'"' + (-l)^a(M/2) if M is even. 



where k does not denote an integer but the wave number 
We now define the matrix 

M-l 



(TM)ifc = ^ V^ij0«fc- (C.6) 



(=0 
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Note that for convenience we have used the complex eigenvectors (1C.3I) . while the ma- 
trix ( 1B.20I) is defined in terms of the real eigenvectors of {A—B)[A+B) and {A+B)[A—B). 
However, these are related by the transformations 

4>k^U(})k and Vfc^t^V'fc 

with the same unitary matrix U . This mapping leaves the right-hand side of equation ( ]C.6I) 
unchanged. Therefore, the two matrices ( 1B.20|) and ( 1C.6|) coincide. 
The matrix (]C.6I) now becomes 

27r(l-l/A./) 

2vr ^ \Ak\ 
For M large enough there exists an integer n < M such that 

Therefore, 

n 

lim A,(M) = g (e^') = ^ (a(j) - ^b{j)) e^^'. 

j=-n 

By taking the limit as M ^ oo of the left-hand side of equation fIC.Tp we obtain (lC.2p . □ 



Appendix D. The Riemann constant K 

In this appendix we will show that the Riemann constant K is given by 



i=2 



2n 



As in PI, let Qi, ... , Q2n-i be the zeros of the theta function 6{u!{z)). Then the function 

e{u;{z)-J2^iQj)-K) 

has the same zeros as 6{lj{z)). Therefore, the quotient of these two functions can be written 
as an Abelian integral of a holomorphic 1-form u: 



z 

V. 



Moreover, all the a-periods of v must vanish. Thus, the right hand side of the above 
equation is in fact a constant C: 
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Therefore, we have 

We will now compute the values of uj{Xi) in the basis ai, . . . , a2n-i, • • • , fc2n-i and show 
that the 2n — 1 points A3, ... , X^n-i are the zeros of 9{uj{z)). We have 



'^j('^2fc+l) 




0< J 


<k<2n-l 


l^j(A2fc+l) 


1 

= -2 + 




0<k<j<2n-l 




= 2^j,k- 


-1, < 


j < k <2n 


UJj{X2k) 


1 

= -2 + 


2^j,k-i, 


1 <k<j <2n. 



If we write u^Xi) as 

cu{Xi) = ^Ni + ^UMi, 

then, from the periodicity (16.61) of the theta function, we have 

e{uj{X,)) = exp i-m {Ni, M,)) e{-uj{X,)). 

Since (A''2j+i, M2j+i) are odd for 1 < i < 2n — 1, we see that 6'(ti;(A2i+i)) = and hence the 
g zeros of 9{lj{z)) are the points A3, ... , Xin-i- Therefore, we have 

2n 

i^=-^^(A2,_l). 

i=2 



Appendix E. The cycle basis (110.181) 



In this appendix we will show that the basis defined in fllO.lSp are canonical. First note 
that, by direct computation, it is easy to check that the intersections between the a-cycles 



are zero 



an-j-i ■ an+i = 0, < j,l < r - 1. 



We will now compute the other intersection numbers by induction. 
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First let us compute the intersection numbers between the tilded basis and the untilded 
basis. We have 

On-fc-1 ■ 0.n+j = —Sk,j (E-2) 

On+fe ■ On+j = <^fc,j (E-4) 

2(-l)'=-^ J+'l<A;; (E.5) 
0, < A; < J - 1. 

r _ / 0, 0<A;<j; ,^ 

an-k^i ■ 6n+i - I 2(-l)^^-^ J + 1 < A;. ^^'^^ 

ttn+k ■ bn-j-1 = (E.7) 

On+fc ■ ^n+j = Sk,j (E-8) 

~ / -1, < A; < J -1; 

bn+k ■ ttn-j-l = 0, j < k. ^ ^ 



On^k~l ■ fln-i-l - ^ 0, J <k. 

, ~ _ / < A; < j; 

bn+k ■ an+j - I 0, J + 1 < A;. 

, ~ _ / < A: < j; 

■ an+j - I 0, J + 1 < A;. 

, r _ f 1, 0<A;<j-l; 

■ bn-,-i - I 0, J < A;. 

, r _ f 1, 0<A;<j; 

o„_fc_i ■ o„_,„i - I J + 1 < A;. 

, r _ r 0, 0< A;<j; 

o„_fc_i ■ On+j - <y J + 1 < A;. 

^n+k ' bn+j = 

where j, k range from to r — 1. 
Now, we have 

Then from (E.9)-(E.16), we obtain the following intersection numbers: 

bn+r-l ■ 5j = 6n+r-l,j, ' bj = 0. 

Next, from (110.181) we have 

bn+k + bn+k-1 = bn+k + + ^^n-fc-l " 26„_fc_i. A; = 1, . . . , r — 1 



E.IO) 

E.ll) 

E.12) 

E.13) 

E.14) 

E.15) 
E.16) 
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From this relation and equation (1E.1I) -(E.16). we obtain 



'n+fc + On+k-l 



~Sj,n+k ~ Sj,n+k-l 



0. j = l,...,2n-l 



Therefore, if we assume that bn+k has the intersection numbers 



bn+k ■ Cij 
bn+k ■ bj 



^j,n+k 

0, j = l,...,2n-l, 



then bn+k~i will have the intersection numbers 



bn+k-l ■ Clj — 
bn+k~l ■ bj = 

Therefore, by induction we see that 



~Sj,n+k-l 

0, J = l,...,2n- 1. l<k 



bn+k ■ Clj 
bn+k ■ bj 



^j,n+k 

0, J = l, 



,2r2-l, k = 0,...,r-l. (E.17) 
We can now compute the intersection numbers of the bn-k-i- We have 

bn-k-l — bn+k = —(in+k + On+fc- = 0, . . . , r — 1 

Therefore, by using flE.ip - flE.17p we obtain 

bn+k ~ bn-k-l 

bn+k - bn-k-l) - bj = 0, j = 1, . . . , 2r2 - 1 
From flE.17p . we see that the intersection numbers for the bn-k-i are indeed given by 

bn—k—l ' 5ij Sj^n—k—1 

bn-k-l ■ bj = 0, j = 1 



~^j,n+k + ^j,n-k-l 



2n-l, A; = 0, ...,r- 1. 



Appendix F. Solvability of the Wiener-Hopf factoriza- 
tion problem 

We now show that the Wiener-Hopf factorization problem (15. 3p is solvable when (3{X) is 
purely imaginary. 

In other words, we have 
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Theorem 9. The following Riemann-Hilbert problem 

T+{z) = Hz)T^{z), \z\ = l 

*(^' - (-/m.) 1a') f^-^) 

whereT^{z) is holomorphic for \z\ < 1 andT_{z) is holomorphic for \z\ > 1 withT_{oo) = 
1 is solvable when /3(A) G iM. 

Proof. We will use the vanishing lemma to proof this theorem. As in [8j, we need to show 
that a certain singular integral operator is a bijection. 

The solvability of the Riemann-Hilbert problem is related to the bijectivity of a singular 
integral operator. Let C be the Cauchy operator 

and let C+, C_ be its limit on the positive and negative side of the real axis 

CM){z) = \unC{f){z±te), z e E. 

Now, define the singular integral operator C$ as in [8]. 

C$(/) = C7+ (F.2) 

Suppose that J — C$ is invertible in i^^(S), and let fi = {I — C$)~^C+(/ — $~^): then the 
function 

f(^)=J + C((J + /x)(J 
is a solution to the Riemann-Hilbert problem (IF.ll) . In fact, we have 
f+{z) = I + C+il - <^-^) + C^fi = I + fi{z) 

f.{z) = f+{z)-I-fi{z) + <^-\z){I + fi{z)) = <l>-\z)f^{z), \z\ = l, 

where the second equation follows from the identity C+ — C_ = /. 

Therefore, in order to show that (IF. ID is solvable when (3{X) G iH, we need to show 
that / — C$ is invertible in i^^(S). 

Using standard analysis (see, e.g'.,[8]), we can show that the operator C$ is Fredholm 
and has index zero. Therefore we only need to show that its kernel is {0}. 

Suppose that the kernel is non-trivial and let (J — C$)/io = 0. Then the function 

foiz)=C[fioiI-<^~')] 
will solve the Riemann-Hilbert problem (IF.ip . but its asymptotic behavior will be 

Tq{z) = 0{z~^), z — >• oo. 
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This means that the function R{z) = Tq (z~^)To(z), where A"^ is the Hermitian conjugate of 
A, is analytic outside the unit circle and behaves like 0{z~'^) at infinity. Thus, by Cauchy's 
theorem, we have 

J^R_{z)dz = 0. 

By making use of the jump conditions, we obtain 

j^R-{z)dz = l^(f^iz))^(fo{z)) dz 

= J^(T^{z))y{z)(To{z))_dz = (F.3) 

From (16. 2p . we see that the eigenvalues of ^{z) are i{X + 1) and i{X — 1). Therefore the 
matrix ^$^(2;) Hermitian and is either positive definite or negative definite for /3(A) G iM.. 

This means that the boundary value of (^o{z)^ on the unit circle is zero. In particular, 

it implies that To{z) = and hence the kernel of the singular integral operator / — (7$ is 
trivial. This concludes the proof of the theorem. □ 
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